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ABSTRACT 

Context. The formation of high mass stars and clusters occurs in giant molecular clouds. Objects in evolved stages of massive star 
formation such as protostars, hot molecular cores, and ultracompact HII regions have been studied in more detail than earlier, colder 
objects. Further progress thus requires the analysis of the time before massive protostellar objects can be probed by their infrared 
emission. With this in mind, the APEX Telescope Large Area Survey of the whole inner Galactic plane at 870 pm (ATLASGAL) has 
been carried out to provide a global view of cold dust and star formation at submillimetre wavelengths. 

Aims. We derive kinematic distances to a large sample of massive cold dust clumps from their measured line velocities. We estimate 
masses and sizes of ATLASGAL sources, for which the kinematic distance ambiguity is resolved. 

Methods. The ATLASGAL sample is divided into groups of sources, which are located close together, mostly within a radius of 2 
pc, and have velocities in a similar range with a median velocity dispersion of ~ 1 km s -1 . We use NTR, N 2 FL, and CS velocities to 
calculate near and far kinematic distances to those groups. 

Results. We obtain 296 groups of ATLASGAL sources in the first quadrant and 393 groups in the fourth quadrant, which are coherent 
in space and velocity. We analyse HI self-absorption and HI absorption to resolve the kinematic distance ambiguity to 689 complexes 
of submm clumps. They are associated with 12 CO emission probing large-scale structure and 13 CO (1-0) line as well as the 870 pm 

dust continuum on a smaller scale. We obtain a scale height of ~ 28 ± 2 pc and displacement below the Galactic midplane of-7 ± 1 

pc. Within distances from 2 to 18 kpc ATLASGAL clumps have a broad range of gas masses with a median of 1050 M 0 as well as a 
wide distribution of radii Their distribution in galactocentric radii is correlated with spiral arms. 

Conclusions. Using a statistically significant ATLASGAL sample we derive a power-law exponent of -2.2 ±0.1 of the clump 
mass function. This is consistent with the slope derived for clusters and with that of the stellar initial mass function. Examining the 
power-law index for different galactocentric distances and various source samples shows that it is independent of environment and 
evolutionary phase. Fitting the mass-size relationship by a power law gives a slope of 1.76 ± 0.01 for cold sources such as IRDCs and 
warm clumps associated with HII regions. 
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1. Introduction 

1.1. High mass star formation 

Current theory of star formation distinguishes between the for¬ 
mation of stars with low and high (> 8 M 0 ) mass. In spite of 
much effort to understand the process that leads to the formation 
of high mass stars, very little is known about their earliest phases 
in contrast to the better understood case of isolated low mass star 
formation. Massive stars release a great amount of energy into 
the interstellar medium through radiation and stellar winds and 
deposit heavy elements in their late evolutionary phases or via 
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supernova explosions, which leads to a chemical enrichment of 
the interstellar medium. They thus influence their environment, 
may trigger the formation of following generations of stars, a nd 
a lso p lay a dominant role in the evolution of galaxies dKennicutti 
[20051) . To advance the understanding of high mass star forma¬ 
tion, many surveys have been conducted so far. Targeted sur¬ 
veys detected strong radio emission from ultracompact HII re¬ 
gions (UCHIIRs), which are formed when the protostar emits 
ultraviolet radiat ion, thus heating and ionizing the remaining 


molecular cloud (Wood & Chu 

rchwell| 198^ Becker et all 199(1 

Kurtz et al.l [1994t (Walsh et al. 

1991 Hoare et al. 20121 k More- 


over, extensive surveys have been made toward samples that 
were selected at far infrared wavelengths (IRAS) to identif y sites 
of high mass stars at the earliest evolutionary phase (IWalsh et al.l 
Il997t ISzvmczak et al.ll200f)t iMolinari et al1l2000l) . In addition, 
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many new massive star forming regions have been revealed in 


surve y s for 6.7 GHz methanol maser emissi on (Szymczak et al 
2002: ISridharan et al.l l2002t ICaswell et al.l 1201 It Green et al 


2009). 


Significant progress in the study of star formation has also 
been achieved by various large-scale infrared continuum sur¬ 
veys. Dust obscuration of the visible light from stars hampers 
the detection of the stellar content of the inner Galaxy, which 
made surveys at longer wavelengths necessary. The Midcourse 
Space Experiment (MSX; [Price et al.ll2001l) observed the whole 
Galactic plane within Galactic latitudes of ±5° between ~ 6 and 
25 pm with a spatial resolution of ~ 18". It provided the data 
necessary to investigate interstellar dust, young stellar objects as 
well as HII regions and Galactic structure. More recently with 
still increasing resolution and sensitivity the Spitzer Galactic 
Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE; 


IBeniamin et al 


__ 2 0031) and th e MIPS Galactic Plane Survey 

(MIPSGAL: ICarev et alJl2009h followed. Emission from 3 to 8 
pm with a spatial resolution of ~ 2" is observed by GLIMPSE 
over a Galactic latitude range of ±1° and longitudes between 
295° and 65°. As a longer wavelength complement MIPSGAL 
surveys 278 deg 2 of the inner Galactic plane at 24 and 70 pm 
with resolutions of 6" and 18". Surveys at longer wave lengths 
are the Bolocam Galactic Plane Survey (BGPS; [Aguirre et al.l 
1201 lh covering -10° < l < 90° and -0.5° < b < 0.5° at 1.1 
mm and the Herschel Infrared Galactic Plane Survey (Hi-GAL; 
iMolinari et al.ll20fd ) reaching a Galactic longitude of ±60° and 
latitude of ±1° from 70 to 500 pm. Another large -sca le project 
is the Red MSX Source (RMS, iLumsden et allI 20 I 3 I) survey, 
which searches for massive young stellar objects (MYSOs) in 
the whole Galaxy. Multi-wavelength observations are conducted 
to distinguish the MYSOs and UCHIIRs from other red sources, 
which allows statistical studies of these stages of high mass star 
formation. 

All these surveys easily find luminous sources in the Galac¬ 
tic plane. However, the disadvantage of these targeted surveys 
is that they are partial, probing only a particular evolutionary 
stage. For example, UCHIIRs detected by their free-free emis¬ 
sion probe recent massive star formation. Moreover, many suf¬ 
fer from incomplete statistics. To overcome these drawbacks, 
the first unbiased submm continuum survey of the whole in¬ 
ner Galactic disk, the APEX Tele scope L arge Area Survey of the 
Galaxy at 870 pm (ATLASGAL: ISchuller et alJl2009t) was con¬ 
ducted. Using the Large APEX Bolometer Camera (LABOCA), 
observations were made in a Galactic longitude range of ±60° 
and latitude of ±1.5° in order to obtain a statistically relevant 
sample of objects associated with high mass star formation at 
various stages as well as to study their distribution and to com¬ 
pare physical properties of these clumps. The definition of a 
clump given by IWilliams et alJ { 2000 1 is adapted in this article. 
We thus consider clumps to be overdense substructures within 
molecular clouds that are coherent in Galactic longitude, lati¬ 
tude, and radial velocity in molecular line maps. C lumps exhibit 
sizes between 0.3 and 3 pc (iBergin & Tafallal2007l) and can form 
whole clusters of stars in high mass star forming regions. Al¬ 
though the ATLASGAL survey is very important to get a global 
view of star formation at submillimetre wavelengths by identi¬ 
fying all massive clumps forming high mass stars in the inner 
Galaxy, a main limitation that applies to all continuum surveys 
is that the distances to the newly found sources are unknown. 
These are needed to determine important parameters such as 
masses and luminosities of the clumps as well as to analyse the 
spiral structure of the Milky Way. To measure distances towards 
a large sample of high mass star forming clumps located at a vast 


range of distances, we need an efficient method, which is the 
kinematic distance estimate (Wilsonl 1972 : Kolnak et al.l 120031 


Roman-Duval et al.ll2009l:lUrauhart et al.ll2014bl). 


Ammonia observations of ATLASGAL sources with kine¬ 
matic distances in the first quadrant of the Galaxy are presented 
in lWienen et alJ (120121) . where the distance ambiguity is not re¬ 
solved. We now complement this data with new NH 3 measure¬ 
ments towards ATLASGAL sources in the fourth quadrant and 
add a new distance analysis towards the combined sample. Sec- 
tion[2]describes molecular line and HI data, which were used to 
derive kinematic distances to ATLASGAL sources. We give de¬ 
tails about the rotation curve and the computation of the errors in 
the distance in Sect. [4] Section0describes the two methods used 
to resolve the kinematic distance ambiguity, HI self-absorption, 
and HI absorption. We show our results in Sect. [ 6 ] and discuss 
them in Sect. [7] A summary of our analysis is presented in Sect. 

m 


1.2. Kinematic distance 

To obtain an estimate for the kinematic distance requires the 
measurement of the radial velocity of a source, which can be 
associated with its galactoce ntric radius using a model of the 
Galactic rotation curve (e.g. [Reid et all 20091 3 ; [Brand & Blitzl 
Il993l) . While the derivation of distances from velocities is 
straigthforward in the outer Galaxy, it is more challenging for 
sources with Galactic radii smaller than that of the Sun, thus lo¬ 
cated in the inner Galaxy. The longitude and velocity of an object 
determine a unique galactocentric radius, but in general allow for 
two different values for its distance to the Sun. Furthermore the 
near and far kinematic distances are equally spaced on the near 
and far side of the point, where the line of sight is tangent to 
the orbit of the object and its radial velocity equals the circular 
velocity. The near and far distances are only at the tangent point 
the same. 

There are many studies, which tried different methods to 
solve the kinematic distance ambiguity (KDA). In addition to 
the derivation of distances to high mass star forming regions, 
pulsar distances are mostly determined using the kinematic dis¬ 
tance method. Measurements of HI absorption at 21 cm against 
bright pulsars leading to lower andjipper distance limits to 
these objects are show n e.g. by iKoribalski et al.l d 19951) and 
[Frail & Weisbergl dl990l) . A common approach used for HII re¬ 
gions consists of measuring an absorption spectrum against the 
continuum free-free emission radiated by the HII region. Near 
and far distances can be distinguished by comparing observed 
absorption line velocities between the source velocity and the 
velocity of the tangent point. Previous work analysed H 2 CO ab¬ 
sorption agdns^_ 6 _^m_coiUinunm emissionto solve the KDA 
dWilsonll 19721 : 1 Watson et alJl2003USewilo et akll2004l) . Since the 
HI abundance is high towards molecular clouds, KDA solutions 
using this probe are more reliable than those resulting from 
H 2 CO absorption. Distances to small samples of HII regions are 
deter mined using HI absorption against r adio continuum emis¬ 
sion dKolnak et alJl2003t iFish et alJl2003h . To re solve the KDA 
to com pact HII regions from the RMS survey, lUrauhart et al. 
(2012) used the same technique. Recent work bv iGarcfa et al. 
(2014) distinguished between near and far distances to groups of 
molecular clumps embedded in the same giant molecular cloud 
with one of them harbouring an HII region using continuum ab¬ 
sorption as well. Because this method requires an HII region to 
be embedded in the molecular cloud, it is only helpful to resolve 
the KDA for a subsample of ATLASGAL sources. However, we 
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also need a technique appropriate for all molecular clumps inde¬ 
pendent of them harbouring an HII region. 

In addition to the HI absorption, we investigate HI self¬ 
absorption (HISA) toward observed NH 3 emission of ATLAS¬ 
GAL sources to resolve the KDA. Using this technique we 
analyse if cold HI in a molecular cloud absorbs warm HI line 
emission in the interstell ar medium. Earli er studies have also 
used this method, e.g. iJackson et alJ ( 2002i t found HISA to¬ 
ward the molecular cloud GRSMC 45.6+0.3, which reveals its 
near kinematic distance. Recently, lUrauhart et al.l (i2014bl) ex¬ 
amined HI absorpti on to determine distances to ~ 800 RMS 
sources. In addition. iRoman-Duval et all d2009l ) analysed HISA 
together with the HI absorption to resolve the KDA toward 
750 molecular clouds observed within the Boston University- 
Five College Radio Astronomy Observatory Galactic Ring Sur¬ 
vey (GRS, [Jackso n et al.l [20061) . Both techniques are also used 
by iBusfield et al. ( 20061) to determine dista nces to a sample of 
massiv e young stellar objects as well as by lAnderson & Banial 
(l2009i) for a sample of HII regions. They investigated the con¬ 
fidence of HISA and both fin d that it agrees with HI a bsorp- 
tion in ~ 80% of their samples. lAnderson & Banial (20091) point 
out several reasons why the HI absorption is still more certain 
for the analysis of distances to HII regions. In addition, reli¬ 
able KDA solutions can be derived from interferometric observa¬ 
tions. The se avo id contamination from large-scale HI emission 
(Urquhart et al.l[ 20 T 2 h and also avoid confusion caused by sev¬ 
eral regions in the beam, although most ATLASGAL sources are 
associated with a single HII region and only 23 are associate d 
with multiple HII regions (see Fig. 3 in lUrauhart et aUl2013bl) . 
However, interferometric observations of a large ATLASGAL 
sample, to which we derive kinematic distances, are prohibitive. 

This paper uses a consistent method to determine kinematic 
distances to a large sample of ATLASGAL sources in the first 
and fourth quadrant. We are aware that our method to derive 
kinematic distances can result in large distance errors for indi¬ 
vidual sources. However, our conclusions do not rely on correct 
individual distances. We focus on the analysis of statistics of pa¬ 
rameters derived from the distances, which is not affected by the 
distance errors, because they average out for our large sample of 
ATLASGAL sources. 


2. Data sets 

2.1. Observations 

We mainly use our subsample of ATLASGAL sources, which 
was observed in the NH 3 (1,1), (2,2), and (3,3) inversion lines 
within a Galactic longitude range from 5° to 60° to derive kine¬ 
matic distances. We measured NH 3 radial velocities of an AT¬ 
LASGAL sample in the first quadrant described in lWienen et al.l 
( 120121 ) and complement it with ammonia observations of a flux- 
limited sample of ATLASGAL sources in the fourth quadrant 
down to about 1.2 Jy/beam (Wienen et al. in prep.). For the mea¬ 
surements in the first quadrant we used the Effelsberg 100-m 
telescope with a beamwidth (FWHM) of 40" at the frequencies 
of the NH 3 (1,1) to (3,3) lines at ~ 24 GHz and a spectral reso¬ 
lution of about 0.7 km s _I . Sources in the fourth quadrant were 
observed using the Parkes 64-m telescope with a 13mm receiver 
as frontend ranging between 16 and 26 GHz and a beamwidth 
of 60" at ~ 24 GHz. We had a Digital Filter Bank (DFB3) as 
spectrometer with a spectral resolution of about 0.4 km s~'. We 
observed two polarizations of the NH 3 (1,1) to (3,3) lines simul¬ 
taneously in position-switching mode. The sourc e sel ectio n as 
well as data reduction are described in I Wienen et al . (2012:1) and 



Galactic Longitude (deg) 


Fig. 1: Galactic longitudes of ATLASGAL sources are compared 
with the radial velo cities w ith CO (1-0) emission in the back¬ 
ground (Dame et al .1120011). The green line indicates the 5 kpc 
molecular ring ( Si mon et al .120060 . 


in Wienen et al. (in prep.). Because we observed only half the 
number of sources in NH 3 in the fourth quadrant than in the first 
quadrant resulting in 315 radial velocities in the fourth quadrant 
compared to 752 velocities in the first quadrant, we added other 
molecular line data for the sources in the fourth quadrant (see 
Sect. l2.2l i. 


2.2. Molecular line data 


A subsample of ATLASGAL sources was followed up in molec¬ 
ular lines. From 1067 ATLASGAL sources observed in NH 3 HI 
data were available for 1056 clumps and we thus used radial ve¬ 
locities of 749 ATLASGAL sources in the first quadrant and of 
307 clumps in the fourth quadrant from NH 3 line measurements 
to derive kinematic distances. We added N 2 H + (1-0) observa¬ 
tions of ATLASGAL sources in the fourth quadrant (Wyrowski 
et al. in prep.). In addition to our own data we use v elocities from 
CS (2-1) measurements of IRDCs from lJackson_gUal] _(j2008(). 
from C S (2-1) observations of UCHIIRs from Bronfman et alJ 
d!996h . and from 13 CO measurements of MYSOs from the RMS 
survey (iUrauhart et al.l(2007 ). We thus obtain radial velocities 
of a total of 1065 sources in the fourth quadrant. The different 
molecular lines are summarized in Table [Q which contains the 
molecular probe, the number of all sources observed in NH 3 , 
and of clumps added from N 2 H + , CS, and l3 CO measurements, 
the number of sources, for which HI data are available and 
which were used to resolve the KDA, their longitude ranges, 
and a reference. The distribution of the clumps used to deter¬ 
mine ki nematic distances is overlaid on the CO (1-0) emis¬ 
sion (Dame et al.l 12001 ) in Fig. Q] The straight line illustrates 
the 5 kp c mo lecular ring, where star formation is actively going 
on (Simon et alil2006l) . Most clumps show strong CO emission, 
which traces the large-scale structure. 


2.3. Archival data 

The VLA Galactic Plane Survey 

To look for HI self-absorption and HI absorption of the ATLAS¬ 
GAL sample in the first quadrant we mainly use the VLA Galac- 
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Table 1: Summary of molecular lines used to obtain radial velocities. 


molecular tracer 

number of 
all sources 

number of 
sources with HI data 

longitude range 
(deg) 

Reference 

NH 3 (1,1) 

752 

749 

5-60 

Wienen et al. (2012) 

NH 3 (1,1) 

315 

307 

300 - 360 

Wienen et al. (in prep.) 

N 2 H + (1 -0) 

301 

301 

300 - 360 

Wyrowski et al. (in prep.) 

CS (2 - 1) 

146 

146 

301 -358 

Jackson et al. (2008) 

CS (2 - 1) 

111 

111 

300 - 355 

Bronfman et al. (1996) 

13 CO(l -0), 13 CO (2-1) 

200 

200 

300 - 350 

Urquhart et al. (2007) 


tic Plane Survey tVGPS. IStil etaDbOObh . The HI 21 cm line 
cubes were created from the survey, which observed the Galactic 
longitude range from 18° to 67° and latitude varying from ±1.3° 
to ±2.6°. The VGPS has an angular resolution of 1', a spectral 
resolution of 1.56 km s -1 , and an rms noise of 2 K per 0.824 
km s -1 channel. 

The Southern Galactic Plane Survey 

For sources with Galactic longitude 1 < 18° HI data are 

availab le from_the Southern Galactic Plane Survey (SGPS, 
iMcClure-Griffiths et al.l2005i) . The whole survey was conducted 
between 1=253° and 358° as well as from 1=5° to 20° with a 
latitude of ±1.5°. The SGPS uses the Australia Telescope Com¬ 
pact Array and the Parkes telescope. The angular resolution of 
the HI data is 2', the spectral resolution 0.8 km s _1 , and the rms 
sensitivity 1.6 K. 


3. Grouping of ATLASGAL sources with known 
velocities 

We describe in this section our method to group sources that 
cluster in position- velocity space according to the friends-of- 
friends algorithm dHuchra & Gelled Il982t IMoore et alJ 1 1993b 
iBerlind et al.112006 ). As a first step we search around each AT¬ 
LASGAL source within a maximum distance and velocity inter¬ 
val for associated sources. We accumulate ATLASGAL sources 
to the same group until no new source that obeys the two criteria 
to the neighbouring sources is found. All distances between at 
least two individual sources within one group must therefore be 
smaller than the maximum distance and the velocities of these 
must be within a given velocity interval. As a next step we vary 
the maximum distance and the velocity interval to analyse the 
effect of these input parameters on the grouping, which leads to 
different divisions of the ATLASGAL sources into groups. More 
details about that are given in Appendix [A] These divisions are 
investigated to reach a better mapping between our grouping and 
properties of known molecular cloud complexes (see Appendix 
ED>. We refer to sources, which are gathered together according 
to their location and kinematics, as a group and denote its spa¬ 
tial extent as a complex, when we compare its structure with the 
ATLASGAL dust continuum emission. An example of such a 
complex within 22.1° < l < 22.6° and 0.25° < b < 0.5° is given 
in the top panel of Fig. [2] The submm dust continuum emission 
is plotted in the background and traces the filamentary structure 
within this molecular cloud complex. It contains three ATLAS¬ 
GAL sources with velocities around 84 km s 1 known from NH 3 , 
N 2 H + or CS observations overlaid as white dots and ATLAS¬ 
GAL sources without known velocities illustrated as green dots. 
The red contoursshow the 13 CO emission from the GRS survey 
(ilackson et al.ll2006h at the velocity of the complex illustrating 


that the ATLASGAL clumps belong to one complex, which is 
coherent in space and velocity. 

Using our method we obtain 296 groups in the first quad¬ 
rant and 393 groups in the fourth quadrant. The group number, 
the mean position, number of sources per group, the mean ve¬ 
locity, the velocity dispersion, and the size of the groups are 
given in Table [2j No velocity dispersion and size are calculated 
for sources, which are not associated with another ATLASGAL 
source within 0.3° and 10 km s _I . 


4. Kinematic distance estimation 

We calculated kinematic distances to ATLASGAL sources, 
for which we had determined the LSR velocity, using the 
iBrand & Blitzl (1 19931) rotation curve, which assumes 8.5 kpc for 
the distance of the Sun to the Galactic centre, Ro, and 220 km s 1 
for the circular rotation velocity at the position of the Sun, © 0 . 

A more recent model using the revised rotatio n parame ters of 
the Milky Way has been presented bv lReid et al.l ( 12009bl) . They 
fit their trigonometric parallax measurements of masers in high 
mass star forming regions to a model of the Galaxy and give 
their results for flat_and more complex rotation curves in Table 
3 of lReid et al l d2009bl) . A comparison of measured HI terminal 
vel ocities i n the fourth quadrant of the Milky Way observed by 
IMcClure-Griffiths & Dickevl (120071) with tangent point velocities 
calculated from the different rotation curves revealed that_the 
measu rements are best fitted by the model from IBrand & Blitzl 
dl993l) . Moreover. iReid et~afl (l2009bl) fit rotation curves to mea¬ 
sured parallaxes of masers in 18 high mass star forming regions, 
which are observed as part of the BeSSeL survey only in the first 
and second quadrants of the Galaxy so far. IReid et akl d2009ai) 
aim at extending their measurements to obtain distances in the 
first and fourth quadrant. Of course by now more sources have 
been observeiQ The ATLASGAL sources are better represented 
by the IBrand & Blitzl d 19931) rotation curve, before trigonomet¬ 
ric parallax measurements will improve the rotation curves in 
the inner Galaxy. To be consistent for all ATLASGAL clumps 
we use one rotation curve for the whole Galaxy. 

We get two solutions from the rotation curve for every ra¬ 
dial velocity of sources inside the solar circle, corresponding 
to a near and far distance. For objects outside the solar circle 
only one distance is physically possible. For sources located at 
the tangent point the near and far distance are similar. We give 
the source name, LSR velocity, the KDA solution (see Sect. [3. 
the kinematic distance calculated from the LSR velocity of each 
source, the distance of the group, in which the clump is located, 
its errors, the gas mass, radius (see Sect. 16. 3L and the group num¬ 
ber in Table [3] The distance of the group is computed from the 
mean of all source coordinates and velocities inside the group. 
To avoid any influence of peculiar motion we recommend the 

1 http://www3.mpifr-bonn.mpg.de/staff/abrunthaler/BeSSeL/index.shtml 
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Fig. 2: a) Example of a complex in the first quadrant: The 870 /rm dust continuum is shown in the background probing the small- 
scale molecular cloud structure. We illustrate ATLASGAL sources with observed velocities of ~ 84 km s _I as white dots and 
ATLASGAL sources without known velocities as green dots. The 13 CO emission at the velocity of the complex is indicated as red 
contours, b) The panel on the left displays the HI emission map of the region shown in a: The HI intensity is integrated in the 
velocity interval around the group velocity and in intervals at smaller and larger velocities. The subtraction of these maps (see Sect. 
153 reveals HI self-absorption plotted in the background with the 870 /;m dust continuum overlaid as blue contours, the l3 CO (1-0) 
emission as green contours, and the 12 CO intensity as white contours. We show ATLASGAL sources with measured velocities as 
red triangles and ATLASGAL sources without velocities as green triangles. We extracted HI lines at the location of observed source 
and at one offset position, which is indicated by the star on the map. The HI ”on-source” spectra are illustrated as solid lines and 
the ”off-source” spectrum as dashed line in the panel on the right. The labels on the top display the galactic coordinates of observed 
source positions. Vertical lines in the spectra indicate the variation of the source velocity resulting from the velocity dispersion of 
clouds or streaming motions. The difference of the on-source and off-source spectra is shown below the HI spectra, c) Example of 
another region, which illustrates HI emission. The HI map in the background is produced and overlaid with the same contours and 
symbols as described for panel b. 
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Table 2: Properties of complexes (see Sect.[3}- The full table is available at CDS. 


Complex 

mean longitude 
(deg) 

mean latitude 
(deg) 

number of sources 
per complex 

mean velocity 
(km s -1 ) 

mean velocity dispersion 
(km s -1 ) 

size 

(pc) 

1 

5.2732 

0.0662 

9 

9.9 

1.4 

22 

2 

5.6161 

-0.0872 

2 

-25.7 

- 

- 

3 

5.8834 

-0.4295 

22 

11.8 

3.9 

28 

4 

5.9901 

-1.2114 

13 

12.3 

2 

104 

5 

6.1909 

-0.3626 

2 

-33.8 

- 

- 

6 

6.2288 

-0.079 

8 

14.8 

3.3 

39 

7 

6.6718 

-0.203 

16 

16.7 

5.3 

22 

8 

7.1655 

0.132 

1 

80.9 

- 

- 

9 

7.3017 

-0.1621 

4 

22 

3.5 

47 

10 

7.2911 

-0.5397 

4 

19.2 

1.5 

3 

11 

7.3395 

-0.01 

2 

31.5 

- 

3 

12 

7.5809 

-0.1494 

5 

153.2 

1.1 

23 

13 

8.0089 

-0.2355 

5 

40 

0.8 

26 

14 

8.1533 

0.2662 

9 

18.7 

0.1 

13 

15 

8.5578 

-0.3362 

17 

37.4 

1.1 

57 

16 

8.4959 

-0.9814 

2 

16 

- 

0.4 

17 

8.9549 

-0.5346 

1 

20.2 

- 

- 

18 

9.0374 

-0.521 

1 

37.1 

- 

- 

19 

9.2132 

-0.1991 

1 

42.1 

- 

- 

20 

9.8551 

-0.732 

4 

27.6 

- 

5 


distance of the group, which is used to determine gas masses and 
sizes. Two of the sources (l,b = 5.62,-0.08;/,/? = 6.19,-0.36) 
have velocities inconsistent with the rotation model and we do 
not assign a distance to them. 

4.1. Calculation of errors in kinematic distances 

Errors in the distance computation originate from the deviation 
of the source velocity from the values one calculates assuming 
circular rotation. This difference is produced by the random ve¬ 
locity dispersion within molecular clouds and in addition sys¬ 
tematic streaming motions, when gas flows through spiral arms. 
We adopt 5 km s 1 given by Brand &Blitzl(l 19 931) as ave rage for 
the ran dom part found in literature. iMcClure-Griffiths & Dickevl 
( 2007 ) determined 6.6 km s 1 for the standard d eviatio n of 
velocity residuals between SGPS data and the iBrand & Blitz! 
(119931) rotation curve to estimate the amount of streaming mo¬ 
tions in the inner Galaxy. Using the random and systematic con¬ 
tributions to non-circ ular motio ns we calculate the error in the 
radial velocity 6\ = V5 2 + 6.6 2 km s _1 = 8.28 km s -1 . We vary 
the LSR velocity of each clump by <5v, the distance uncertainty 
is then the difference between the distance caculated with vlsr 
and the distance obtained with vlsr ± 8.28 km s _1 . To deter¬ 
mine errors in distances to the groups we take into account that 
the random velocity dispersion decreases for groups with many 
sources and vary the mean radial velocity of a group, < vlsr >, 
by 


<Wg r o U p = J (6.6 km s" 1 ) 2 + (cry / Viv) (1) 

with the velocity dispersion of sources in a group, cr v , and the 
number of clumps in a group, N. For the calculation of the er¬ 
ror we distinguish between ATLASGAL sources, which are lo¬ 
cated close to the tangent point, and those that are not. In the 
first quadrant ~ 10% of the groups (~ 5% in the fourth quad¬ 
rant) have a mean radial velocity slightly smaller than the tan¬ 
gent point velocity, but the velocity used for the error calcula¬ 


tion (< vlsr >+ ^ v group) is larger than the tangent point ve¬ 
locity. Approximately 5% of the groups in the first and fourth 
quadrant have mean radial velocities slightly larger than the tan¬ 
gent point velocity, which likely results from streaming motions, 
while the velocity used to estimate the distance uncertainty (< 
V LSR > —hvgroup) is still smaller than the tangent point velocity. 
Approximately 1 % of the groups in the first and fourth quadrant 
have mean radial velocities and < vlsr > -<Wg roU p greater than 
the tangent point velocity. These ATLASGAL sources, which lie 
close to the tangent point, have a near and far distance similar to 
the distance computed with the tangent point velocity. 

Most sources, 85% in the first quadrant and 89% in the fourth 
quadrant, are not located close to the tangent point. The near dis¬ 
tance can vary by a large amount for clumps with radial veloc¬ 
ities close to 0 km s 1 as a result of streaming motions. Hence, 
the near distance is unknown for the subsample with velocities 
between ~ -6 km s -1 and ~ 8 km s' 1 in the first quadrant and 

from-9 km s 1 to ~ 6 km s -1 in the fourth quadrant. There are 

a few sources with extreme velocities: Some clumps in the fourth 
quadrant with longitude at ~ 355° have radial velocities around 
100 km s -1 . These are located in the surrounding of Bania’s 
Clump 1 (iBaniall9771) . which have mostly a non-circulm'veloc- 
ity and for which we do not determine a distance. ICaswell et al.l 
(l2010t) consider that they might be in the Galactic bar because of 
their high radial velocities. 

Further estimates of distanc es to HII regions in the extreme 
inner Galaxy were made bv lJones et al.l (120131) . They use HI ab¬ 
sorption to probe the longitude-velocity distribution of features 
within a Galactic longitude range up to 10° and give lower limits 
to the line-of-sight distance by comparing the HI absorption ve¬ 
locity with the radio recombination line velocity of HII regions. 

5. Resolving the Kinematic Distance Ambiguity 

Two methods, one using HI self-absorption and the other HI ab¬ 
sorption, gives us an indication whether to chose the near or the 
far distance of ATLASGAL sources. While HISA can be used 
for all molecular clumps, the HI absorption needs a background 
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Table 3: Kinematic distances to ATLASGAL sources determined using HISA and HI absorption, gas masses, and radii (see Sect. 
Q). The full table is available at CDS. 


Name 

v(l,l) 

KDA 

Distance 

Distance comp i ex 

Mgas 

radius 

Complex 


(km s- 1 ) 

solution 

(kpc) 

(kpc) 

(Mo) 

(pc) 


G5.35+0.10 

11.21 

n 

3.04 

0 '7C+1.07 

Z,,/ -1.67 

2.72 

0.31 

1 

G5.39+0.19 

10.58 

n 

2.91 

0 '7C+1.07 
-1.67 

2.35 

0.19 

1 

G5.64+0.24 

7.94 

n 

2.28 

0 '7C+1.07 

Z,,/ -1.67 

3.01 

0.81 

1 

G5.62-0.08 

-25.74 

f 

- 

- 

- 

- 

2 

G5.83-0.51 

16.22 

n 

3.64 

9 Q4+0.96 

-1.42 

3.24 

1 

3 

G5.83-0.40 

7.64 

n 

2.16 

9 Q4+0.96 
-1.42 

2.80 

0.78 

3 

G5.89-0.39 

9.18 

n 

2.47 

9 Q4+0.96 
-1.42 

3.14 

0.64 

3 

G5.89-0.29 

10.11 

n 

2.65 

9 Q4+0.96 
-1.42 

2.74 

0.56 

3 

G5.89-0.32 

10.11 

n 

2.65 

9 94+0.96 
-1.42 

2.65 

0.58 

3 

G5.90-0.44 

9.1 

n 

2.45 

9 94+0.96 
-1.42 

2.99 

0.71 

3 

G5.90-0.43 

6.64 

n 

1.91 

9 94+0.96 
-1.42 

3.11 

0.66 

3 

G5.91-0.54 

14.79 

n 

3.41 

9 94+0.96 
-1.42 

2.46 

0.36 

3 

G6.13-0.63 

15.85 

n 

3.49 

O Q4+0-96 
-1.42 

2.94 

1 

3 

G6.21-0.59 

18.38 

n 

3.78 

9 94+0.96 
-1.42 

- 

- 

3 


HII region to be associated with the ATLASGAL source. We 
resolve the kinematic distance ambiguity by combining our re¬ 
sults from the two techniques, which are described in the next 
sections. Similar studies have already analysed the two meth¬ 
ods to estimate distances to smaller samples or mostly limited 
to the first quadrant of t he Galaxy dAnderson & Banial l2009t 
iRoman-Duval et akll2009h . Our KDA resolutions to a large sam¬ 
ple of ATLASGAL sources in the first and fourth quadrant give 
an unbiased 3D view of massive star forming clumps in the 
Galaxy. The high number of kinematic distances allows us to 
derive important parameters such as masses and sizes of molec¬ 
ular clouds in the first and fourth quadrant and to study their 
statistics. 

Based on the groups found in Sect. [3] we define complexes 
on their association with extended 12 CO and 13 CO emission and 
association with ATLASGAL sources without known velocities. 


5.1. HI self-absorption 

The HI gas within dense molecular clouds is cold with a tem¬ 
perature of about 10 K because it is shielded from the external 
interstellar radiation field. In contrast, warm atomic hydrogen is 
dis tribu ted throughout the ISM and has a temperature of ~ 100 
K (IClark et al.1119620 . When emission from the warm HI back¬ 
ground at the far distance is absorbed by cold foreground HI 
at the same velocity, HI self-absorption occurs. This results in 
an absorption line in the HI 21 cm spectrum at the velocity of 
the molecular tracers observed towards the ATLASGAL clumps 
such as NH 3 , N 2 H + or CS. A strong absorption therefore hints 
at the near kinematic distance. Because the diffuse warm HI gas 
is abundant in the whole Galactic plane, it is also located in the 
background of a cloud at the far distance. However, the radial 
velocity increases with the distance only up to the tangent point 
and decreases afterwards. The warm HI background is therefore 
at lower velocities than the clump, which leads to no HI self¬ 
absorption toward this source. In addition, foreground radiation 


is emitted by warm HI gas at the same velocity as the observed 
clump. If HISA is not seen, the probability of the far kinematic 
distance is therefore increased. 

However, this method is affected by a few uncertainties. 
It relies on the approximation of a uniform distribution of HI 
gas, although an absorption might also arise from a lack of HI 
emissio n in some spec ific location. Previous HISA studies (e.g. 
[Baker & Burtonlll979l) found criteria, which increase the confi¬ 
dence of the technique, such as the limitation of the Full Width 
at Half Maximum of the absorption to a few km s -1 (< 5 km s -1 ) 
or absorption lines should have an intensity of at le ast 10 K to 
separa te HISA from background variations. lAnderson & Banial 
(2009) translated this temperature requirement into a column 
density of 6 x 10 21 H 2 molecules needed to have the minimum 
HI column density (cm -2 ) for HISA. The ATLASGAL clumps 
fulfill this_criterion becaus e for the sources that are observed in 
NfL. iWienen et all (i2012l) obtain values of the H 2 column den¬ 
sity between 2.6 x 10 21 and 1.8 x 10 23 cirT 2 with an average of 
1.7 X 10 22 cm -2 , which lies above this column density limit. 

5.2. HI Self-absorption analysis 

To examine if HISA exists toward a complex of ATLASGAL 
sources, we create maps of average HI intensity at the veloc¬ 
ity of the complex and look for HI absorption features toward 
source positions with known velocities. If absorption is found, 
this indicates the near kinematic distance of sources within the 
complex. The absence of any absorption favours the far kine¬ 
matic distance solution. To increase the HI intensity contrast in 
case of an absorption, we first average the HI emission in the ve¬ 
locity range of sources in a complex (the ”on” map), from which 
we subtract an HI map containing the HI intensity in a velocity 
interval different from that of the complex (the ’’off” map). To 
create the HI on map, we average the HI 21 cm line cubes over 
the velocity range v± Av with the mean radial velocity of sources 
in the complex, v, and the mean of their linewidths, Av. Two off 
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maps are built by averaging the HI intensity over two times the 
mean linewidth at a sligthly smaller and larger velocity than the 
mean radial velocity of the complex, over the velocity intervals 
v — Av — 4 ± Av km s _l and v +Av + 4 + Av km s _1 . The two maps 
at the offset velocities are averaged and then subtracted from the 
on map. 

In addition, we extract the HI 21 cm spectrum for source po¬ 
sitions in each complex that show a NH 3 , NTH ! or CS emission 
peak (on spectrum), and compare it to the HI 21 cm line at an 
offset position in the same complex, free from emission of the 
mentioned molecular tracers at the velocity of the complex. One 
would reveal HISA as an absorption line in the HI on spectrum 
in contrast to the spectrum at the offset location. However, an 
absorption might also result from a fluctuation of the HI back¬ 
ground at the source velocity, which adds some uncertainty. We 
therefore use this spectral analysis of HISA as additional infor¬ 
mation to the HI average intensity map, the HI spectrum is ad¬ 
vantageous if it is not clear from the HI map whether it shows an 
absorption or not. If the KDA resolutions of both disagree, we 
decide for the distance assigned by the analysis of the HI map. 

Panel b) of Fig. [2] shows the HI map of the complex plotted 
in panel a), it results from the subtraction of the maps at offset 
velocities from the map at the mean radial velocity of the com¬ 
plex. The clumps with observed velocities around 84 km s _1 are 
indicated as red triangles and sources without known velocities 
as green triangles. To define the boundary of a c omp lex we over¬ 
laid 12 CO (1 - 0) emission from lDame et alJ (12001 4 tracing the 
large-scale structure on the HI map, illustrated as white contour 
in Fig. [2] The 13 CO (1-0) line and the 870 pm dust continuum 
are shown as green and blue contours to probe molecular clumps, 
which are identified by ATLASGAL sources, on a smaller scale. 
We assign ATLASGAL sources without known velocities, which 
are traced by contours of 12 CO emission, to the same complex. 
We give the mean of the observed velocities of sources within a 
group as determined in Sect. [3] to the sources, which have orig¬ 
inally no measured velocity. We plot HI spectra of the region 
shown in panel b) on the right as solid lines, which are extracted 
at the position of three observed sources. The HI spectrum at 
one offset position, indicated by the star on the map, is added as 
dashed line to the HI spectra. The coordinates of the extracted 
HI lines are given above the spectra. The vertical lines around 
the source velocity of the HI spectra indicate its variation due to 
velocity dispersion of clouds and streaming motions. The differ¬ 
ence of the on-source and off-source spectra is illustrated below 
the HI lines. Because absorption is revealed in the map as well 
as in the spectra at the source velocity, the complex in panel b) 
is assigned to the near kinematic distance. 

In addition, another HI map is plotted for a complex within 
30.5° < Z < 31.3° and -0.05° < b < 0.45° in panel c) of Fig. 0 
It contains observed ATLASGAL sources with velocities around 
40 km s~‘ illustrated as red triangles and a large number of AT¬ 
LASGAL sources without known velocities as green triangles. 
We assign the mean velocity of the complex of 39.81 km s 1 
(see Table [2} to the sources with no measured velocity. The HI 
map of this complex reveals HI emission and no absorption is 
visible in the HI spectra as well, which hints at the far kinematic 
distance. 

5.3. HI absorption 

Another technique to resolve the kinematic distance ambiguity is 
the analysis of absorption toward sources with strong radio con¬ 
tinuum emission, e.g. HII regions. Their thermal free-free con¬ 
tinuum emission at 21 cm, which has a higher brightness tem¬ 


perature than that of neutral HI gas in clouds, is absorbed by 
HI in molecular clouds located between the HII region and the 
observer. The absolute value of the radial velocity increases with 
the distance from the Sun up to the tangent point velocity and de¬ 
creases at larger distances. Foreground clouds, which absorb the 
21 cm continuum of an HII region at the near kinematic distance, 
thus have HI absorption lines at velocities up to that of the HII 
region, which is known from observations of NH 3 , NiH + or CS 
used in this article (see Sect. 12.2b in the same molecular cloud. 
For a 21 cm continuum source located at the far kinematic dis¬ 
tance, HI absorption at velocities up to the tangent point velocity 
is revealed. 

IWilsonl d 1972b analysed first H 2 CO absorption against 21 cm 
continuum emission as a tool to resolve the KDA for 28 HII re¬ 
gions. However, H 2 CO is not as broadly distributed as HI and 
probes rather molecular than atomic gas. Since the inner Galaxy 
is filled by a much larger amount with atomic than molecular gas, 
there might be not enough H 2 CO in foreground molecular clouds 
leading to its non-detection near the tangent point. The absorp¬ 
tion spectrum toward an HII region may consequently not exhibit 
high velocities and the near distance might be assumed, although 
the object is located at the far distance. Since the HI abundance 
is high in molecular clouds, KDA sol utions using this probe ar e 
more reliable than the HISA method (iAnderson & Baniall2009l) . 

5.4. HI absorption method 

To utilize this technique we use HI 21 cm line cubes from the 
VGPS and SGPS data sets. To recognize if an absorption dip 
in the HI spectrum is produced by HI absorption toward an HII 
region we calculate the difference between the HI spectrum ex¬ 
tracted at the position of the source and in a region around it. 
Assuming that the HI emission is relatively smoothly distributed 
and extended we subtract the HI line emission and the resulting 
spectrum reveals absorption lines in case of an HII region. To 
derive the HI spectrum in an envelope around the source we av¬ 
erage the emission within 120". We also tried out a larger size of 
420" to extract more HI emission inside a larger region around 
the source. However, this has a slightly different shape than the 
emission coming from the source itself producing HI intensity 
fluctuations. These might be considered as absorption dips with 
velocities exceeding the source velocity and thus indicating the 
far distance instead of the near distance. To avoid those HI inten¬ 
sity fluctuations resulting from the emission in the larger region 
around the source, we decided to use the small envelope. Exam¬ 
ination of the HI spectra reveals that we can only use them to 
resolve the KDA if they satisfy the following criteria: 

1. We estimate the noise level <x from part of the spectrum free 
from absorption and use only spectra with absorption fea¬ 
tures over 5cr. Moreover, different shapes of HI lines at the 
source position and in the surrounding result in intensity fluc¬ 
tuations around 0, when the two spectra are subtracted. A 
dip in the HI spectrum must fall below —1.5 times the max¬ 
imum positive intensity to be recognized as absorption line. 
The fluctuations should not exceed 5<x for sources in the first 
quadrant. This is not set in the fourth quadrant because in¬ 
spection of HI spectra for a sample of known HII regions 
shows that too many sources would not be considered as HII 
regions. 

2. To assure that the source is indeed associated with a strong 
HII region the continuum flux is extracted from 21 cm con¬ 
tinuum emission maps. A source is identified as HII region 
if the mean continuum level is higher than 3 to 5<x, where 
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(a) 

G8.14+0.23 G337.92-0.46 




(b) 


G37.37—0.23 



G337.12-0.18 



(C) 

G49.37 —0.30 



Velocity (km/s) 


G309.89 + 0.40 



Velocity (km/s) 


Fig. 3: Example of HI spectra towards HII regions in the first and fourth quadrant extracted from 21 cm line cubes from the VGPS 
and SGPS. The source velocity is illustrated as solid red vertical line, the tangent point velocity as solid blue vertical line, and the 
velocity, at which the absorption is below 5cr, as dashed-dotted green vertical line. This limit of -5cr is indicated by the solid black 
horizontal line. In addition, the maximum positive intensity of HII regions in the first quadrant must be lower than 5cr, which is 
shown as dashed black horizontal line. The dashed black vertical lines +8.28 km s around the source velocity display its variation 
because of the velocity dispersion of clouds and streaming motions (see Sect. ED. Examples of sources in the first quadrant are 
shown on the left and of clumps in the fourth quadrant on the right, a) There is a large difference between the tangent point velocity 
and the velocity, up to which HI absorption is seen and which is close to the source velocity for G8.14+0.23 and G337.92-0.46, 
which are thus located at the near distance. They are located in the diagonal right region in Fig.[4](red points) confined by the amount 
of streaming motions, b) Our analysis reveals HI absorption at velocities close to the tangent point velocity for G37.37-0.23 and 
G337.12-0.18, which are assigned to the far distance and indicated as blue points in the horizontal right part in Fig. |4] c) We locate 
G49.37-0.30 and G309.89+0.40 at the tangent point because the source velocity is within 8.28 km s -1 around the tangent point 
velocity. They are displayed as purple points in Fig. [4] 


<x is the continuum noise (~ 100 mJy/beam for the sample 
in the first quadrant and ~ 600 mJy/beam for the sources in 
the fourth quadrant). The different flux limits result from dif¬ 
ferent sensitivities and beam sizes of the VGPS and SGPS 
described in Sect. 12.31 


3. There must be an HI absorption line at the velocity of the 
molecular cloud, which harbours the ATLASGAL clump 
consisting of cold gas observed in NH 3 , N 2 H/ or CS. This 
assures that the HII region is also embedded in the molecular 
cloud. 
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We obtain 65 HII regions in the first quadrant and 146 HII re¬ 
gions in the fourth quadrant, for which we can distinguish be¬ 
tween near and far distances using the method from lKolnak et al.l 
( 2003 1). They determined the source velocity (v s ), the velocity of 
the tangent point (v t ), and the maximum velocity of absorption 
lines in the 21 cm continuum (v a ) of a sample of HII regions 
located in the first quadrant. Their comparison of v t - v s and 
v t - v a showed that the sources are mostly located in two re¬ 
gions, where v t - v a is around 0 and where v t - v a is increasing 
with ri_sing Vj-Vj. To assign kinematic distances to those regions, 
iKolpak et all (120031) ran simulations with 10000 HII regions to 
analyse the velocity differences, which yields that sources with 
v t - v a close to 0 are located at the far kinematic distance, while 
those with growing v t - v a and v t - v s are at the near distance. 

We obtain the source velocity of the HII regions from NH 3 , 
N 2 H + , and CS emission associated with the dust clumps (see 
Sect. 12.21) . The tangent point velocity is calculated using the rota¬ 
tion curve bv lBrand & Blitz! (119931) for the longitude range from 
- 20 ° to 60°, while it is compute d in the fourth quadrant with 
the linear function by IMcClure-Griffiths & Picked d2007l) . who 
fitted HI terminal velocities as a function of Galactic longitude 
from 300° to ~ 339°. Figure Q] shows HI absorption line spectra 
of ATLASGAL sources in the first and fourth quadrant placed 
at the near and far kinematic distance and at the tangent point. 
The tangent point velocity is indicated as solid blue vertical line 
and the source velocity as solid red vertical line. We determine 
the noise level of the spectra, cr, and give the limit of 5cr, which 
should not be exceeded by the maximum positive intensity of the 
sources in the first quadrant, as dashed black horizontal line. We 
measure the maximum velocity in the first quadrant at which the 
absorption is below 5cr, which is illustrated by the solid black 
horizontal line, and lower than —1.5 times the maximum pos¬ 
itive intensity. The same is considered in the fourth quadrant, 
where the absorption velocity is the minimum velocity, which is 
indicated by the dashed-dotted green vertical line in Fig. [3 The 
velocity dispersion of clouds and streaming motions, which re¬ 
sult in differences of the source velocity from circular rotation, 
are shown as dashed black vertical lin es ± 8.28 km s_|_ around 
the source velocity (lBurtonll~l 97 if, IMcClure-Griffiths & Dickevl 
120071) . Seventeen sources of our sample in the first quadrant 
and seven clumps in the fourth quadrant have a source velocity 
within 8.28 km s ~ 1 of the tangent point velocity or lie beyond the 
tangent point, for which we assign the distance of the tangent. 

To distinguish between near and far distance we plot the ve¬ 
locity differences v t - v a against v t - v s for our subsamples of 
HII regions in the first and fourth quadrant in Fig. [4] as shown 
in IKolnak et al.l (12003}) for simulated HII regions. They studied 
the quality of their KDA resolutions, which shows that 90% of 
sources in the diagonal right part bounded by solid red lines in 
Fig. |4] are at the near distance and 90% of HII regions in the 
horizontal right region, confined by solid red lines as well, are 
at the far distance. Those sources are assigned to distances with 
a high degree of confidence. The dashed red lines around them 
indicate additional space taking int o accou nt streaming motions 
of 8.28 km s _1 . According to IKolnak et al.l (120031) the triangle 
at the lower left part of the plot contains sources, for which the 
resolution of the kinematic distance ambiguity is more uncer¬ 
tain. The dotted line in the triangle separates 21 cm continuum 
sources at the near distance in the upper region from those at the 
far distance below the dotted line with a lower degree of confi¬ 
dence. The velocity differences of identified HII regions in the 
first quadrant are shown as black points in Fig. [4] and those of 
HII regions in the fourth quadrant in green, for which we show 
the inverted velocity differences to be able to present the two 



Fig. 4: Comparison of the differences between the tangent point 
velocity (v t ), the maximum velocity of absorption lines in the 
21 cm continuum (v a ), and the source velocity (v s ) of HII re¬ 
gions. The subsample in the first quadrant is shown in black and 
sources in the fourth quadrant in green. To be able to display 
the two in one plot, we invert the velocity differences of HII re¬ 
gions in the fourth quadrant. The solid and dotted lines confine 
regions, where 21 cm continuum sources are assig ned near and 
far distances using the method from IKolnak et al.l (120031) . The 
examples of ATLASGAL sources shown in Fig.[3]are indicated: 
Clumps at the near distance are shown in red, sources at the far 
distance in blue, and clumps at the tangent point in purple. 


samples in one plot. Six sources in the first quadrant and 34 
sources in the fourth quadrant lie between the horizontal and 
diagonal regions. In the sample from iUrauhart et aD (120121) 10 
HII regions are also in this location in the same plot of velocity 
differences. They took their distance estimate from the litera¬ 
ture, mostly from spectrophotometric measurements and H 2 CO 
absorption, which suggests the near distance for these sources. 
As a consequence we adopt their KDA resolution and place the 
sources from our sample between the horizontal and diagonal 
regions in Fig. [4]at the near distance. 


5.5. Methods to resolve the KDA 


To obtain distances we first apply the two methods, HISA and 
HI absorption, separately and then combine the results. This is 
illustrated as flow chart in Fig. [3 We summarize here the most 
important steps. 


Before we look for HISA we divide the ATLASGAL sam¬ 
ple into groups of sources, which are located close together and 
have similar velocities (see Sect. [3]). Then, we derive maps of av¬ 
erage HI intensity to investigate HISA toward the groups. These 
are overlaid with molecular line emission as contours tracing 
small- and large-scale cloud structure, which allows us to assign 
ATLASGAL sources without known velocities to the molecu¬ 
lar cloud complexes and thus to determine distances to them as 
well. In addition, we look for absorption lines in the HI spec¬ 
trum at source positions with known velocities. This is the first 
step of the flow chart in Fig. [3 together with the investigation 
of absorption lines in the 21 cm continuum extracted as de¬ 
scribed in Sect . 15.41 As next step in Fig. |3 we use the method 
from lKolnak et al. ( 20031 ) for sources showing an absorption fea¬ 
ture, explained in Sect. 15.41 to distinguish between near and far 
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distances. For clumps not associated with HI absorption we use 
only HISA. 

To combine the two methods we distinguish between com¬ 
plexes that contain one or several continuum sources. The third 
step in the flow chart shows the case of only one continuum 
source in a complex. If the KDA resolutions from HISA and HI 
absorption agree, we will use the assigned distance, indicated as 
dashed red arrows in Fig. [5]for the near distance and solid green 
arrows for the far distance. If the far distance is determined from 
HISA and the near distance from HI absorption as shown by the 
solid black arrows, the HI intensity map might not show self¬ 
absorption because star formation can be ongoing in the cloud 
and emission of protostars heats their environment, which can be 
detected at 8 pm by the GLIMPSE survey. To identify these cases 
we look for 8 pm emission in the surroundings of our source and 
will assign the near distance if we find GLIMPSE counterparts 
for sources with a high S/N ratio of the absorption line in the 
21 cm continuum. If there are no associations with GLIMPSE 
sources, the far distance is chosen. For the near distance deter¬ 
mined from HISA and the far distance from HI absorption as 
shown by the dashed blue arrows, we use the far distance be¬ 
cause the self-absorption might not result from absorption of the 
HI background, but from the absorption of the HII region em¬ 
bedded in our source. 

For complexes that have more than one 21 cm continuum 
source one can distinguish between near and far distance if the 
KDA resolutions from absorption features in the 21 cm contin¬ 
uum agree. If they do not agree and e.g. two HI spectra give 
a near and far distance from HI absorption, while HISA deter¬ 
mines the far distance, one has to examine the S/N ratio of the HI 
spectra. If one of those is favoured because of a higher S/N ratio, 
one can resolve the KDA, otherwise no distance is assigned. 

5.6. Resolved kinematic distances 

As described in Sect. l5.5l we first arrange our sample with known 
velocities of 749 sources in the first quadrant and 1065 clumps 
in the fourth quadrant into groups that are coherent in space and 
velocity (see Table[I}. Comparison of these complexes with line 
and dust continuum emission reveals molecular cloud structure, 
which is important to associate ATLASGAL sources without 
known velocities with the complexes as well. Using HISA we 
distinguish between near and far distances for the whole AT¬ 
LASGAL sample including 1489 sources in the first quadrant 
and 2069 clumps in the fourth quadrant with and without known 
velocities. We arrange 165 ATLASGAL sources in the first quad¬ 
rant and 184 clumps in the fourth quadrant in more than one 
complex and without velocity information we cannot determine 
distances to them. 

For sources identified as HII regions with strong continuum 
emission we obtain additional indication for near and far dis¬ 
tances from our analy sis o f absorpti on in the 21 cm continuum 
following the work bv lKolnak et al.l ( 2003 1. While we locate all 
sources with |v t - v s [ < 8.28 km s -1 at the tangent point taki ng 
non-circular motions into account (see Sect. 15.41) . IKolnak et al.l 
(2003) place only HII regions with v t - v s < 0 at the tangent 
point. Based on their criterium we find 36 sources with a 21 cm 
continuum in the first quadrant from our sample with known ve¬ 
locities at the near distance, 23 at the far distance, and 6 at the 
tangent point, 107 HII regions in the fourth quadrant are at the 
near distance, 38 at the far distance, and 1 at the tangent point. 
Considering non-circular motions the number of sources at the 
tangent point is increased to 17 in the first quadrant, where we 
obtain 31 near and 17 far HII regions. We assign a near distance 


to 110 sources with a 21 cm continuum in the fourth quadrant, a 
far distance to 29, and 7 are placed at the tangent point. 

A comparison of near and far distances to HII regions de¬ 
termined from HISA as well as HI absorption using the method 
from iKolpak et al.l (120031) yields an agreement in the results of 
the two methods of 58% in the first quadrant and 71% in the 
fourth quadrant. The relatively low agreement is expected be¬ 
cause HISA gives better results for the near distance and HI ab¬ 
sorption for the far distance (see Sect. 1531 ). 

We combine our KDA resolutions obtained from HISA and 
HI absorption as described in Sect. 15.51 Sources with velocities 
within 8.28 km s -1 of the tangent point velocity are placed at the 
tangent point. The distance assignments of clumps with mea¬ 
sured velocities in the first and fourth quadrant are summarized 
in Table [4] These ATLASGAL sources are divided into com¬ 
plexes, which we locate at the tangent point if the difference of 
the mean radial velocity and the tangent point velocity is smaller 
than the amount of non-circular motions estimated by equation 
Q] The resulting distances to complexes are also shown in Ta¬ 
ble [4] A subsample of 27 ATLASGAL sources with resolved 
kine matic distances has also been observed by the BeSSeL sur¬ 
vey (iReid et al.ll20 14h to derive trigonometric parallax distances. 
These are compared with the kinematic distances in Appendix 
□ 

Using flux limited ATLASGAL subsamples in the first and 
fourth quadrant, which goes down to a peak brightness of 1.5 
Jy/beam, we obtain a ratio of near-to-far kinematic distances 
of 1.5 ± 0.21 in the first quadrant and of 2.4 + 0.35 in the 
fourth quadrant. The distributions in the first and fourth quad¬ 
rant agree within a 3<x error. T his trend is similar to the result 
from lRoman-Duval et al.l d2009l) , who obtain the near distance to 
2.3 times more molecular clouds revealed within the GRS survey 
than the far distances. The KDA r esol ution for samples of HII re- 
gions from lUrauhart et all (12012l) and lAnderson & Banial (120091) 
yields a larger number of sources located at the far distance than 
at the near distance. Since the submm continuum and 13 CO emis¬ 
sion are trac ing all molecular clumps, o ur source sample as well 
as that from lRoman-Duval et al.l (120091) c ontain a larger number 
of n earby, low mass clumps. I n contrast, llJrauhart et al.l (12012l) 
and lAnderson & Ba nia (2009) are probing preferably massive 
star forming clumps that are typically found at larger distances. 

Recently, lEllsworth-Bowers et al.l ( 20131) determined a dis¬ 
tance probability density function (DPDF) from a kinematic 
distance as the likelihood and prior probabilities related to an¬ 
cillary data sets to resolve the KDA to BGPS_clumps. This 
work is expanded bv lEllsworth-Bowers et al.l (120141) . who intro¬ 
duce two new prior DBDFs relating BGPS clumps to objects 
with trigo nometric parallax distances from the BeSSeL survey 
(iReid et al.ll2009b1) and to HII regions with resolved kinematic 
distances using HI absor ptio n methods from the HII Region Dis¬ 
covery Surveys (HRDS. lBania et akll20ldl2012 ). A comparison 
of the kinematic distanc es to ATLA SGAL sources with the dis¬ 
tances presented in lEllsworth-Bowers et al. (20141) reveals that 
the agreement is ~ 70% within 2 kpc. The remaining fraction 
of inconsistent distances do not affect our statistical analysis be¬ 
cause distance errors average out for our large ATLASGAL sam¬ 
ple. 


6. Results 

With the resolution of the kinematic distance ambiguity of a 
large sample of 3558 ATLASGAL sources we study their Galac¬ 
tic distribution and analyse the large-scale structure of the Milky 
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Fig. 5: Flow chart illustrating HI self-absorption and HI absorption to resolve the KDA. The upper panel shows the main points of 
the two methods and the lower panel combines their KDA resolutions. 


Table 4: Summary of the KDA resolution of ATLASGAL sources and complexes. 


Sample 

number of sources at 
near distance 

far distance 

tangent 

point 

Total 

number of complexes at 
near distance 

far distance 

tangent 

point 

Total 

first quadrant 

383 

332 

34 

749 

143 

148 

5 

296 

fourth quadrant 

737 

290 

38 

1065 

234 

141 

18 

393 


Way. In addition, we use the distances to calculate important 
physical properties such as masses and sizes. 

6.1. Galactocentric distribution of ATLASGAL sources in the 
Milky Way 

We show the galactocentric distribution of the ATLASGAL 
clumps in the first and fourth quadrant in the upper and lower 
panel of Fig. [6] 

6.1.1. Galactocentric distribution in the first quadrant 

The top panel indicates two peaks at a galactocentric radius, 
RoaU between 4 and 5 kpc and at ~ 6 kpc. A comparison of 
these peaks with the large-scale structure of the Milky Way us- 
ing a model with four spiral arms (e.g. iGeorgelin & Georgelinl 
119761 ) reveals that the enhancement of sources at R(- ia | = 4.5 kpc 
can be attributed to the intersection of the Scutum-Centaurus arm 
and the Galactic bar. The largest fraction of sources contributing 
to that peak, 19%, are assigned to complex 172 (l = 30.767°, b = 
-0.050°), which contains the star forming region W43 (see Ap- 
ncndix lfTTl i. The enhancement at 6 kpc coincides with the galac¬ 
tocentric radius of the Sagittarius arm. Most sources from that 
peak are assigned to two complexes, ~ 27% are located in com¬ 
plex 50 (/ = 14.833°, b = -0.602°), the high mass star forming 
region M16, and 19% are in complex 261 (Z = 49.2°, b = -0.7°), 
which is the massive giant molecular cloud W51 (see Appendix 

IB3 . 

The two peaks in the Galactocentric dist ributio n of the first 
quadra nt are also detected in HII regions by lAnderson & Banial 
(2009). in Bolocam Galactic Plane survey (BG PS) sour ces 
obs erved in the 1.1 mm dust continuum (iDunham et al.l 
201J, and in a samp le of 6.7-GHz methanol masers by 
Green & McClure-Griffithsl (1201 ll) . However, an additional 


peak at a galactocentric radius of ~ 8 kpc is also obtained 
by previous s t udies, e.g. in a sam ple of RMS s ources by 
lUrauhart et al.l ( 201 1). in IRDCs by iJackson et al.l d2008l) . in 
high mass star forming regions with FIR colours of UCHIIRs 
dBronfman et al.l 120001) , and in mo lecular clouds observed in 
13 CO by the GRS survey dRathborne et alJ 120091) . This peak 
results from the buildin g up of sources on the solar circle, which 
are primarily nearby. iPaladini et al.l d2004l) mention that the 
high concentration of HII regions detected at R( ja | = 8 kpc is 
removed in the galactocentric distribution of sources with LSR 
velocities larger than 110| km s -1 . The comparison of the number 
distribution of RMS_sources with radial velocities greater than 
|12| km s _1 (iUrquhart et al.l I 20 TTI) to that of the whole sample 
also reveals a lack of the peak at ~ 8 kpc. The small number of 
ATLASGAL sources in the first quadrant at this galactocentric 
radius might consequently result from a few observations of 
clumps with small velocities. 


6.1.2. Galactocentric distribution in the fourth quadrant 

The histogram including the whole ATLASGAL sample in the 
fourth quadrant reveals three peaks at galactocentric radii of ~ 3 
kpc, between 5 and 6 kpc and at approximately 7 kpc. The en¬ 
hancement of sources at 3 kpc can be attributed to the near 3-kpc 
arm and the southern end of the Galactic bar. Most clumps at 3 
kpc, ~ 20%, are in complex 363 (/ = 351.653 °,b = 0.176°), 
where IRAS 17200-3550 is located. It was observed at 5 GHz 
with the VLA by iBecker et al.l (11994 and is likely a UCHIIR. 
The enhancement between 5 and 6 kpc corresponds to the Norma 
arm and Scutum-Centaurus arm. A large subsample, 31% of 
the number of clumps between 5 and 6 kpc, clustering at the 
Scutum-Centaurus arm is in complex 210 (/ = 333.19 °,b = 
-0.36°), which contains the G333 giant molecular clouds (see 
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Appendix IB.21 ). The peak at ~ 7 kpc can be attributed to the 
Scutum-Centaurus arm as well. Most sources, 42%, are assigned 'T' 
to complex 377 (/ = 353.115°, b = 0.720°), which harbours the $ 
HII region G353.186+0.8 l ocate d at the near kinematic distance 8 
from HI absorption by lOuireza et al.l (120061) . This is in agree- S 
ment with our KDA resolution using HI self-absorption and HI 's 
absorption. Another large contribution to the number of clumps, j> 
32%, originates from complex 356 (Z = 351.317°, b = 0.661°), J 
which is the NGC 6334 molecular cloud complex. Some of “ 
the peaks in the galactocentric distribution of the ATLASGAL g 
clumps in the fourth quadrant are similar to features obtained by Q 
previous studies. A samp le of 6.7- GHz methanol masers in the 8 
fourth quadrant ( Green & McClure-Griffithsl201 ll) also revealed ^ 
the two peaks at 3 - 3.5 kpc and 6 - 6.5 kpc. 

6.1.3. Differences in the Galactocentric distribution between 
the subsamples in the first and fourth quadrant 

The RMS survey did not observe the region within |/| < 10° 
and is thus not sensitive to sources within 4 kpc. Consequently, 
lUrquhart et al.l (120121) obtain only one strong peak between 5.5 
and 6 kpc for a sample of compact HII regions from the RMS 
survey in the fourth quadrant in contr ast to three peaks for their 
sample in the first quadrant ( lUrquhart et al.ll 20 Tll) . This led them 
to deduce that HII regions are not located symmetrically in the 
first and fourth quadrant. Taking the ATLASGAL sample in the 
fourth quadrant into account we detect three peaks, while there 
are two concentrations of sources in the first quadrant at ~ 4 and 
6 kpc. The comparison of the histograms for the subsamples in 
the first and fourth quadrant therefore indicates a difference in 
the overall galactocentric distributions as well, which is consis¬ 
tent with spiral arm structure rather than a molecular ring. 

6.2. Scale height and distance from the Galactic midplane 

We calculate the distance from the Galactic plane in pc from 
the Galactic latitude at the distance of the clumps and show the 
number distribution of the whole ATLASGAL sample in the first 
and fourth quadrant as a black histogram in Fig. [7] An exponen¬ 
tial fit to the data gives a scale height of 27.93 ± 1.97 pc. We fit 
the height distributions of the subsamples in the first and fourth 
quadrant as well, which results in scale heights of 23.62 ± 1.05 
pc in the first quadrant and 30.63 ± 3.32 pc in the fourth quad¬ 
rant. The comparison of the two samples shows that their scale 
heights are consistent within a 3 cr uncertainty. The red histogram 
in Fig. [7] gives the number distribution with Galactic latitude of 
ATLASGAL sources with masses above the completeness limit 
of 1000 M 0 (see Sect. 16. 3I >. A fit to those data yields a scale 
height of 32.99 ± 2.6 pc, similar to that of the whole ATLAS¬ 
GAL sample. The scale height derived from the latitude dis¬ 
tribution of high mass star forming clumps is therefore not af¬ 
fected by their masses. The comparison of our scale height to 
other samples in the inner Galaxy reveals that we obtain a sim¬ 
ilar scale height_a_s derived from 6.7-GHz methanol masers of 
27 ± 1 pc bv lGreen & McClure-Griffithsl (1201 ll) and ~ 30 pc de¬ 
termined fro m a sampl e from the Arecibo Methanol Maser Sur¬ 
vey (iPandian et al.l2009h . The scale height of the whole ATLAS¬ 
GAL sample also agrees with that determined for extended radio 
sources in the first quadrant, which are likely compact and ultra¬ 
compa ct HII regions, of 25 pc dBecker et al.ll99(jl) . lBeuther et al.l 
(12012 1) estimate an upper limit of the scale height of 46 pc for 
ATLASGAL submm clumps, which is higher than our measure, 
because they do not have distances and only approximate them 




galactocentric distance (kpc) 

Fig. 6: Number distribution of ATLASGAL sources with galac¬ 
tocentric radii, the clumps in the first quadrant are plotted at the 
top and the sources in the fourth quadrant at the bottom. 


using a sample of YSOs. A similar scale height of 47 pc is deter¬ 
mined for an ATLASGAL sample of embedded, compact clumps 
by ICsengeri et al.l d2014l) . This was calculated without distance 
information of the individual sources and therefore resulted in 
a larger scale height than our estimate. The scale height of AT¬ 
LASGAL clumps above 1000 M Q (see Sect. l6.3l > agrees with that 
derived fro m a R M S sam ple in the first quadrant of 30.2 ± 2.8 
pc dUrauhart et al.li20Tlh . A fit to the latitude distribution of the 
whole RMS samp le reveals a larger scale height of 37.7 ± 0.8 
pc ( llJrauhart et al.ll2014bh because it also includes sources in the 
outer Galaxy, which can be located at larger distances from the 
Galactic plane. lUrquhart et al.l d2013bl) obtain a scale height of 
20.7 ± 1.7 pc for compact and UCHIIRs, which lies slightly be¬ 
low our value estimated for the whole ATLASGAL sampl e. A 
larger scale height of 29.1 + 3 pc is derived by lUrquhart et al.l 
( 2013b! ' for the submm host clumps, in which the compact and 
UCHIIRs are still deeply embedded, which agrees with our value 
measured for the whole ATLASGAL sample. 

In addition, we obtain the distance of ATLASGAL sources to 
the Galactic midplane from the exponential fit to the histograms 
illustrated in Fig. [7] The peak of the distribution of the whole AT¬ 
LASGAL sample is below the Galactic midplane at -6.68 ±0.97 
pc. A similar displacement of-10.34 + 0.51 pc and -4.09+ 1.66 
pc is also derived for the subsamples in the first and fourth quad¬ 
rant. This is in agreement with the peak of the latitude dis- 
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tribution of ATLASGAL sources in the longitude range from 
-30° to 20 °, which is skewed to a negative value (ISchuller et al.l 
120091) . Our results are consistent with previous studies, which 
also locate their samples below the Gala ctic mi dplane us- 
ing e.g. 1.1 mm clumps from the BGPS (iRosolowskv et al.l 
I 2010 I) . an AT LAS GAL sampl e consisting of embedded, com¬ 
pact sources dCsengeri et al.l20f4l). bubbles in the inner Galactic 
disk identified with GLIMP SE (IChurchwell et al.ll2006l) . UCHI- 
IRs dBronfman et al. 20001). and CO observations o f molecu¬ 
lar clouds ( Sanders et al.l 19841 : IClemens et al.lll988l) . It is usu¬ 
ally assumed that this offset resu lts from the location of the 
Sun above the Gala ctic plane (e.g. iHumphrevs & Larsenll 1 9951 
ISchuller et alll2009l) . 



Fig. 7: Number distribution of ATLASGAL subsamples in the 
first and fourth quadrant with the height above Galactic plane is 
plotted in black. The histogram of sources with masses above the 
completeness limit of 1000 M e is shown in red. 


6.3. Mass and size distribution 


The gas mass of the ATLASGAL cl umps is derived from the 
dust emission dKauffmann et al.l2008l) 


M gas (M G ) = 1.2 x 10 


-14 


S870/imd 3 d 2 ^exp^ 144 *}° 4 j - lj 

~idj) 


( 2 ) 


with t he 87 0 pm flux dens i ty, Sxi o inn , given in IContreras et alJ 
(l2013h and llJrauhart et al.l (1201 4al) : the wavelength, A, in pm; 
the distance to the clump, d, in kpc; and the dust temperature, Tj, 
with Tf\ = 7’| <nl under the assumption of equal gas and dust tem¬ 
peratures. We average the kinetic temperatures obtained from 
NH 3 observations and compute masses with one temperature of 
20.8+2.9 K for the subsample in the first quadrant and one value 
of 23.1+5.5 K for the sources in the fourth quadrant. An ab¬ 
sorption coefficient, ,_k, of 1.85 cm 2 Ip, at 870 ^m at a gas density 
n(H) = 10 6 cur 3 (lOssenkonf & Hennindl 19941) is used resulting 
from models of dust grains with thick ice mantles and the ra¬ 
tio of metallicity to the solar metallicity, Z/Z Q , is assum ed to be 
L These parameters are equivalent to those from ISchuller et all 
(120091) . Un certainties in the gas mass estimates are discussed in 
Sect. 5.2 in lwienen et afl (120121) . 

The gas masses are plotted against kinematic distances for 
the ATLASGAL sources in the first quadrant as black triangles 
and the clumps in the fourth quadrant as red points in the lower 
panel of Fig. [ 8 ] The error bars indicate the average distance 


and mass errors. We also added a contour plot by dividing the 
kinematic distance range into bins of 2 kpc and the logarithm 
of the mass range into bins of 0.25 and counted the number of 
sources in each bin (see upper panel of Fig. [8j. The broad dis¬ 
tribution of kinematic distances with a peak at 4 kpc and gas 
masses with a peak at 630 M e is similar for the two subsamples. 
Sources with small masses exhibit a velocity close to 0 km s _1 
and are therefore located at near distances. However, the vari¬ 
ation of the velocity about 8.28 km s” 1 (see Sect. 14.11 ) results 
in large errors in the distance and in the mass of these sources. 
Fig. [8] reveals a trend of increasing mass with distance result¬ 
ing in part from an enhanced number of clumps, which may fall 
within the beam at larger distances. The curve shown in Fig. [8] 
indicates the mass corresponding to a minimum 5<x detection for 
point sources. Within the whole distance range the effective com¬ 
pleteness of the survey to compact unresolved sources is ~ 600 
M 0 . Approximately 65% of the ATLASGAL sources in the first 
and fourth quadrant with derived gas masses are above 600 M G , 
which is ap pro ximat ely the limit for a clump to form high mass 
stars dCsengeri et al.l20l4l) . However, the majority of sources are 
extended with respect to the beam, their masses are significantly 
above the catalogue point source sensitivity limit and therefore 
the turn over in the mass distrib ution, which is ~ 1000 M e (see 
Fig. |9) (lUrquhart et al.ll2013al) . is a better estimate of the true 
limit of the mass sensitivity of ATLASGAL. 

The rotation curve bv lBrand & Blitz : ( 1993 !) yields large dis¬ 
tances of ~ 20 kpc for one source with l = 348.551°, b = 
-0.339° and a velocity of 10.4 km s 1 and 24 kpc for another 
clump with l = 349.721°, b = 0.121° and a velocity of 17.6 
km s -1 . Because these velocities also allow us to locate the 
sources in the Far 3 kpc arm, we give a more realistic distance 
of 11 kpc. 

Some extreme sources in Fig. [8] are e.g. the two clumps 
(/ = 338.374°, b = -0.152° and l = 338.406°, b = -0.205°) lo¬ 
cated far away at a distance of 16.47 kpc and with gas masses of 
8.86xl0 3 and 2.12 x 10 3 M 0 . They are in the same complex, for 
which we can use HI self-absorption as well as HI absorption to¬ 
ward the more massive clump associated with an HII region. We 
obtain the far distance, which is_in agreement with the kinematic 
distance given by iRusseill ( 2003 ). Details about other sources, 
which have extremely high masses, are discussed in Appendix 

ED 

To estimate the size distribution of the ATLASGAL sample 
we take the source radius, which is no t deconvolved from the 
beam, from IContreras et al.l (1201 3l) and lUrquhart et al.l (l2014al) 
and converted it into pc with the distance of the source. The dis¬ 
tribution of radii of the clumps in the first and fourth quadrant is 
similar, ranging from 0.01 to 2.7 pc with a peak at ~ 0.3 pc. They 
are compared with kinematic distances for the ATLASGAL sam¬ 
ple in the first quadrant as black triangles and the sources in 
the fourth quadrant as red points in Fig. [TO] The upper curve 
indicates radii calculated usin g 1.3', abov e which large-scale 
emission is filtered out (ISchuller et all I2009I) . while the lower 
curve shows radii, which correspond to the APEX beam radius 
of 9" at the wavelength of 870 pm. ATLASGAL sources are 
extended relative to the beam and we obtain a trend of increas¬ 
ing source radius with rising distance. As already revealed in the 
mass distribution. Fig. [TO] shows that the ATLASGAL sample 
basically consists of mole cular clumps, which ex hibit radii be¬ 
tween ~ 0.15 and 1.5 pc (iBergin & Tafallal[2007l) , and clusters 
containing several clumps. 
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Fig. 8 : Correlation plot of the kinematic distances and gas 
masses for the ATLASGAL sources in the first quadrant as black 
triangles and the clumps in the fourth quadrant as red points. 
The curve shows a minimum 5cr detection. The upper panel il¬ 
lustrates the contour plot, for which we counted the number of 
sources in each kinematic distance bin of 2 kpc and each log¬ 
arithmic gas mass bin of 0.25. The contours give 10 to 90% in 
steps of 20% of the peak source number per bin. The error bar 
displays the average distance and mass errors. 



Fig. 9: Histogram of the logarithm of the gas mass. The vertical 
line displays the turnover in the mass distribution. 


Fig. 10: Dependence of the source radius not deconvolved from 
the beam on the kinematic distance to the ATLASGAL sample in 
the first quadrant as black triangles and the clumps in the fourth 
quadrant as red points is shown as scatter plot in the lower panel. 
The lower curve indicates radii corresponding to the beam radius 
of 9", the upper curve illustrates radii computed with 1.3'. The 
binning of the kinematic distance in the contour plot in the upper 
panel is 1.5 kpc and that of the logarithm of the radius is 0.5. The 
contour levels from 10 to 90% are shown in intervals of 20% 
of the peak source number per bin. The error bar indicates the 
average distance and radius errors. 


6.4. Clump mass function calculation 

With the gas masses calculated in Sect. 16.31 we can examine the 
statistics of clump masses. We first analyse statistically the influ¬ 
ence of the assumption of a constant kinetic temperature, used 
in the mass calculation, on the slope of the clump mass function. 
The distribution of gas masses resulting from the observed ki¬ 
netic temperature of each source is thus compared with masses 
computed with the mean kinetic temperature of the sample. To 
investigate if the two are statistically different, we perform a 
Kolmogorov-Smirnov (KS) test using the ATLASGAL subsam¬ 
ple in the first quadrant, which was observed in NH 3 . The KS 
test does not contradict that the two clump mass distributions are 
the same. Cumulative distribution plots in Fig.|TT]reveal similar 
mass distributions, which are calculated using observed kinetic 
temperatures and a constant temperature. To derive the clump 
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mass function we can therefore use gas masses computed with 
an average kinetic temperature (see Sect. [ 63 - 

Two approaches are usually used to characterize the clump 
mass distribution, which are the differential and cumulative 
clump mass function. The upper panel of Fig. IT 2 l illustrates the 
differential mass distribution of sources with distances between 
2 and 5 kpc, which is fitted by a powe r law as a function of mass 
bins having equal logarithmic widths dKennicutt & Evani[2012l) 

dN/d\og\o(M) = AM y (3) 

with the number of clumps, dN, in a mass bin of width d\ogio(M) 
and constants A and y. The power-law exponent y is related 
to the index a in the power law fitting the number distribution 
against mass, which is of ten used to describe structures in gas, 
dKennicutt & Evansil2012|) 

dN/dM = AM a (4) 

by a — y — 1. We divide the mass range of the whole ATLAS- 
GAL sample at any distance into different numbers of bins and fit 
differential clump mass functions to the mass range above 1000 
M g , for which we are complete (see Sect. 16.31) . The best fit to 
the data gives similar values for a ranging from -1.75 to -1.89 
for varying bin widths with an average of -1.82. We also deter¬ 
mined the_slope using a method based on maximum likelihood 
dClauset et al.il2009l) . which is independent of the histogram bin¬ 
ning and results in a = -1.82 + 0.02 consistent within the er¬ 
rors with the steeper values obtained by our analysis. We also 
used different bin numbers to plot the mass distribution of AT- 
LASGAL sources with heliocentric distances from 2 to 5 kpc, 
which is selected around the peak in the heliocentric distance 
histogram. As illustrated in Fig. [ 8 ] we probe a large range in 
masses over a small distance range. Fitting the differential clump 
mass function above 1000 M Q , which indicates the turnover for 
sources within 2 to 5 kpc, leads again to a narrow range of a be¬ 
tween -2.19 and -2.38 for different bin widths with an average 
of -2.29, which is similar to the values for a obtained for the 
whole ATLASGAL sample at any distance. The method based 
on maximum likelihood (IClauset et al.1120091) gives a slope of 
-2.26 ± 0.05, which is in agreement with the power-law expo¬ 
nents from our investigation. The upper panel of Fig.[12]gives an 
example of the differential clump mass function fitted to masses 
of the sources from 2 to 5 kpc resulting in a power-law exponent 
of -2.19 + 0.11. 

We also fit the mass distribution of the ATLASGAL subsam¬ 
ple between 2 a nd 5 kpc by a lognormal function, which is given 
by equation 5 in lPeretto & Fulled (1201 Oh 

diotlm = Aexp (“ ( loglo(M) - lo g‘° K-k)) 2 /2^ 2 ) (5) 

with the constant A, the peak mass L/p ea k, and the disper¬ 
sion, cr. All masses are given in units of M 0 , hence 

is in (logioMo) -1 . The best fit to our data reveals A — 
188 (logio Mq)- 1 , Mpeak = 700 M e , cr = 0.49 (see lower panel 
of Fig. m, and is a better expression of the mass distribution 
than the power law. 

The cumulative mass function describes the fraction of 
clumps with masses greater than a specific value: it is shown 
in Fig. [13] The distribution hints at a steep power law for the 
clumps with high masses and a flatter power law for sources 
with lower masses. We use a double power law for the cumu- 
lativ e mass function with a break mass, of 3000 M 0 

iReid & Wilsonll2006h 

n (> m )=^:; e h :f“ iow v» |o w- (6) 
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where a\ ow and ffhigh are the power-law exponents below and 
above M b rea k- Equation [6] is used for masses smaller than the 
break mass and 

N(> M ) = AM" hi § h+1 (7) 

for masses equal to and greater than the break mass. Our fit 
to the mass distribution of sources within 2 to 5 kpc results in 
«low = -1.83 and (Thigh = -3.28. The power-law exponent of 
-3.28, which describes the high mass portion of the mass spec¬ 
trum, is steeper than the values of a, obtained from the fit of the 
differential clump mass function. 

The cumulative mass function of ATLASGAL clumps with 
distances between 2 and 5 kpc is better described by a lognor¬ 
mal distribution than by the double power law. A relation for the 
lognor mal function is given by equation 10 in IReid & Wilsonl 
(20061: 

N(>M) = i(l-erf((logio(M)-A 0 )/V2A 1 )). (8) 

Its fit to our data is illustrated in Fig.[T3]and results in Ao = 6.44 
andAi = 1.15. 

To summarize, our fit of the differential clump mass function 
to the whole ATLASGAL sample as well as to sources within 
2 to 5 kpc leads to consistent power-law exponents, while high 
masses have a steeper index of the cumulative mass function. 
The masses of the ATLASGAL clumps can be fitted best by log¬ 
normal distributions. 



Fig. 11: Cumulative distribution function of the gas mass, which 
is computed using the observed kinetic temperature of each 
source shown as solid black line and using the average kinetic 
temperature of the sample illustrated as dotted red line. 


6.5. Volume density 

The FL volume density is calculated from the relation n = 
3 M gas /(47rR^ ff ) with the gas mass (see Se_ct. 16.3b . M gas, and 
the effective radiu s, R,.(u as given in IContreras et all 620 1 3l) 
and lUrauhart et al.l (l2014al) . A correction of the mean molec¬ 
ular wei ght per Hi molecule for helium abundance of mg, = 
2.72mjj <;A 1 lenl 1 197 3[) is applied. The ATLASGAL clumps ex¬ 
hibit volume densities between 226 cm~ 3 and 3.4 x 10 6 cm~ 3 
with a median of 7842 cm~ 3 . The number distribution of the 
volume dens ity is shown in Fig. [14] Our values are similar to 
those of IPunham et al. (201 i ll, who calculate volume densities 
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Fig. 12: Differential clump mass function of ATLASGAL 
sources with distances from 2 to 5 kpc fitted to masses above 
1000 M e in the upper panel yields a power-law exponent of 
-2.19 + 0.11. The vertical line displays the completeness limit 
of 1000 M 0 . A lognormal function is fitted to the mass distri¬ 
bution of this ATLASGAL subsample resulting in A — 188, 
f^peak = 700 M e , and cr = 0.49. Mp ea k is below the complete¬ 
ness limit and the peak itself is thus not well constrained. 




Fig. 13: A double power law fit to the cumulative mass func¬ 
tion of the subsample between 2 and 5 kpc gives an exponent of 
-1.83 for masses below 3000 M G and of-3.28 for masses above 
that in the upper panel. A lognormal mass distribution is fitted in 
the lower panel, which leads to Aq= 6.44 and A i = 1.15. 


between 120 cm ~ 3 and 8 x 10 5 cm ~ 3 but obtain a lower me¬ 
dian of 1072 cm ~ 3 from 1.1 mm continuum observations towards 
BGPS sources. We are tracing a smaller range of volume densi¬ 
ties than the d ust c lumps observed in the 1.2 mm continuum by 
iBeltran et alJ ( 2006 ). who derive values between 10 3 and 3 x 10 8 
cm ~ 3 with a median of 4 x 10 4 cm~ 3 . In addition, our volume 
densities are larger than those of giant molecular clouds in th e 
first and fourth quadrant (iPame et al.l 1 986i: iGarcfa et al.ll2014l) . 
which lie between 5 and ~ 300 cm~ 3 , because objects identified 
by the ATLASGAL survey are more compact than molecular 
clouds. 

The ATLASGAL sample consists of high mass clumps 
in different evolutionary stages from cold cores in an early 
phase such as IRDCs to more evolved objects with higher 
temperatures, e.g. UCHIIRs, as revealed by a comparison 
of their NH 3 ( 1 , 1 ) linewidths and rotational temperatures 
(IWienen et al.i 120121) . To investigate any differences in the 
volume density distribution we divide the ATLASGAL sources 
in two subsamples: Precluster clumps, which have a faint or 
no mid-infrared or far-infrared counterpart, do not harbour 
an embedded protostar or protocluster, while ATLASGAL 


clumps with embedded or associated luminous infrared sources 
emit at 21 /jm. Because the ratio of the 21 fi m flux from the 
MS X survey to the 870 yum peak flux density as shown in 
IContreras et al.l (1201 3l) estimates the heated part of a source, we 
use this measure to distinguish between cold sources in an early 
phase or warm clumps in a later phase. Co mpariso n with more 
sensitive and higher resolution WISE data dWright et al.l[2010b 
is presented in lCsengeri et al.l (120141) . A histogram of the 21 jjm 
flux to the 870 yum peak flux density ratio reveals a peak at 0.8, 
which divides the ATLASGAL sources into a cold sample with 
a ratio smaller than 0.8 and a warm sample with a ratio greater 
than 0.8 corresponding to UCHIIs and YSOs. A comparison 
of the volume density distribution of the two samples indicates 
no difference between the two with similar median values of 
8334 cm -3 for the cold clumps and of 7753 cm " 3 for the warm 
sources. This is in agreement with the analyis of [Beltran et al.l 
( 2006 ). who derived similar FL volume density distributions 
of 1.2 mm dust clumps showing MSX emission and sources 
without MSX counterpart. 
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Fig. 14: Histogram of the H 2 volume density of the ATLASGAL 
sample. Most sources exhibit a volume density of 7842 cm -3 . 


7. Discussion 

7.1. Identification of complexes with known giant molecular 
clouds 

We searched in the literature for the complexes with the largest 
numbers of ATLASGAL sources with measured velocities (see 
Sect. 12.11 and 12.2b given in Table [3 Those had already been 
identified as giant molecular cloud complexes in previous stud¬ 
ies, from which we obtain additional information. Analysis 
of the source distr ibuti on as a function of Galactic longi¬ 
tude dCsengeri et al.ll2014l) already reveals that giant molecular 
clouds exist as statistically significant concentrations of compact 
sources. We summarize our associations in Table ED More de¬ 
tails can be found in the Appendix. 

7.2. Galactic distribution of scale height and molecular gas 
midplane 

We use the ATLASGAL sample in the first and fourth quad¬ 
rant to determine the radial distribution of molecular clumps. We 
plot the distance from the Galactic plane against galactocentric 
radii of the ATLASGAL subsamples in Fig. [15] clumps in the 
first quadrant are illustrated in the top panel and sources in the 
fourth quadrant in the lower panel. Individual complexes with 
extremely high distances from the Galactic mid-plane are indi¬ 
cated as red triangles and discussed in the Appendix lB.il 

We divide the range in galactocentric radius from 3 to 8 kpc 
of the ATLASGAL subsamples in the first and fourth quadrant 
with masses larger than the completeness limit of 1000 M G (see 
Sect. 16.31) as well as the galactocentric radii of the two combined 
into bins of width 1 kpc and plot the number distribution of the 
whole ATLASGAL sample with the height above/below Galac¬ 
tic plane for every bin in galactocentric radius in Fig. [16] An ex¬ 
ponential function is fitted to the distributions, which yields the 
scale height and the distance from the mid-plane in each bin with 
the error given by the rms value. We show the scale height dis¬ 
tribution against galactocentric radius, taken as the mean of each 
bin, in the upper panel of Fig. [17] where the subsample in the 
first quadrant is presented in red, the sources in the fourth quad¬ 
rant in green and the two combined in black. The scale height 
of the clumps in the first quadrant is approximately constant be¬ 
tween 21 ± 3 pc and 33 ± 5 pc from a galactocentric radius of 
3 to 7 kpc and rising to ~ 48 ± 13 pc at 8 kpc. The sample in 
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Fig. 15: Distance from the Galactic mid-plane compared with the 
galactocentric radius for ATLAGAL sources in the first quadrant 
in the top panel and for clumps in the fourth quadrant in the 
bottom panel. Sources with a large displacement from the mid¬ 
plane are indicated as red triangles. 


the fourth quadrant exhibits a scale height distribution varying 
between ~ 8 ± 1 pc and ~ 29 ± 3 pc at 6 kpc and increasing up 
to 44 + 7 pc at larger galactocentric radii. 

The whole ATLASGAL sample shows a slightly rising scale 
height distribution from ~ 18 ± 2 pc to ~ 45 ± 6 pc between 
3 and 8 kpc with approximately constant scale heights from 4.5 
to 6.5 kpc, which are mainly determined by a few complexes 
with large numbers of sources. These are discussed in the Ap- 
nendix lBTl The variation of the scale height with galactocentric 
radius has already been studied by previous work. Observatio ns 
of the CO (1-0) line in giant molecular clouds dSanders et al.l 
11984 1 also show a trend of a slightly increasing scale height be¬ 
tween Real ~ 2.5 and 10 kpc. The scale height distribution of the 
subsamples in the first and fourth quadrant are not significantly 
different, but reveal an approximately constant scale height from 
3 to 8 kpc in agreement with the scale height resulting from_CO 
surveys in the first and fourth quadrant (iBronfman et ak 119881) . 
Using HI data in the inner Galaxy fLockmani d 19841) also obtain 
a narrow range in the scale height compared to the flaring of 
the Galactic midplane i n the outer Galaxy at g alactocentric radii 
larger than 9 to 10 kpc dWouterlo ot et all l990l). Thi s is also con¬ 
sistent with the analysis of Urquhart et al.l d2014bl) . who inves¬ 
tigated 13 CO, CS, and NH 3 lines of a sample of RMS sources. 
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Real = 5 - 6 kpc 


Real = 4 - 5 kpc 


Real = 3 - 4 kpc 
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Fig. 16: Histograms of the height above Galactic plane for every 
bin in galactocentric radius of the subsamples in the first and 
fourth quadrant. The distributions are fitted by an exponential. 



Their scale height is also slightly rising from about 20 pc to ~ 30 
pc between 4 and 8 kpc, while the increase is larger, up to 200 
pc, at a galactocentric radius of 11 kpc. 


7.2.1. Variation of the distance from the mid-plane within the 
inner Galaxy 

The lower panel of Fig. [XT] shows the displacement from the 
mid-plane derived from the exponential fit to the distributions 
plotted in Fig. [16] against the galactocentric radius of ATLAS¬ 
GAL clumps. The clumps in the first and fourth quadrant have 
a similar distribution of the distance from the midplane within 3 
kpc to 8 kpc. Figure [17]reveals that the ATLASGAL sample lies 
close to the Galactic plane, with a distance varying between -20 
pc and 20 pc. This is a flat distribution from 3 kpc to 8 kpc com¬ 
pared with the in crea sing mid-plane displacement of the RMS 
sample ( llJrquhart et a l. 2014b) and of OB stars associated with 
UCHIIRs (iBronfman et alJ2000h up to 200 pc between 9 kpc and 
14 kpc. CO lines in giant molecular clouds (ISanders et al.ll 19841) 
also show an approximately constant distribution of the displace¬ 
ment around -20 pc, although it is restricted to n egative val- 
ues. M olecular clouds observed in CO surveys (IBronfman et al.l 
1988 1 in the first and fourth quadrant are bounded to the Galactic 
plane as well, slightly varying between -40 pc and 60 pc, similar 
to distances of the ATLASGAL sources. 

Moroever, we analyse if there is a dependence of the gas 
mass on the distance from the plane. The height above/below 
Galactic plane of ATLASGAL sources with masses above the 
completeness limit of 1000 M Q is divided into bins of 20 pc 
and the gas mass is averaged within these bins. The distribu¬ 
tion of the mean mass in the first and fourth quadrant is ap¬ 
proximately constant over the range of distances from the mid¬ 
plane. Some clumps with high masses are located at large height 
above or below the Galactic plane. These are in complexes, for 
which we could only use the HI self-absorption method to re¬ 
solve the KDA, because they do not contain any HII region. The 
distance assignment is therefore uncertain and as a consequence 
their masses and displacement from the Galactic plane as well. 
The most massive clumps are expected to be located close to the 
Galactic plane because of the Galactic gravitational potential. 
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Fig. 17: Dependence of scale height on the galactocentric radius 
in the top panel and the distance from the disk of the Milky Way 
shown against galactocentric radii in the lower panel. The sub¬ 
sample in the first quadrant is indicated in red, sources in the 
fourth quadrant in green, and the two combined in black. 


The largest masses are indeed distributed at small distances, but 
many objects of lower masses between 1000 M e and 10 4 M e are 
also located near the mid-plane, which lowers the mean value of 
the gas mass. 


7.3. Mass distribution of ATLASGAL sources in the Milky 
Way 

With resolved kinematic distances we are able to study the 3D 
distribution of a large sample of ATLASGAL clumps. We plot 
the logarithm of the galactocentric radius of sources above the 
completeness limit of 1000 M Q against the azimuth, 8, in Fig. [18] 
Their location is analysed with respect to the large-scale struc¬ 
ture of the Milky Way illustrated by the model by Valleel d 19951) . 
From a statistical examination of various spiral arm models us¬ 
ing different methods such as HI emission, HII regions, magnetic 
fields, CO emission, thermal electron gas IValleel ( 1995 1 built a 
logarithmic four-arm model of the Milky Way with a mean pitch 
angle of 12.5 °. Using this value we illustrate spiral arms given 
by the IValleel d 1995i) model as coloured straight lines in Fig. [18] 
They appear twice on the map because they wrap around, and 
ATLASGAL sources within 2 kpc of the Galactic centre are ex- 
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eluded where the rotation curve is unreliable. The left panel indi¬ 
cates that ATLASGAL clumps are located close to spiral arms as 
well as in the interarm region, the source density is high around 
0 = 0° with clumps at the near distance. The right panel of Fig. 
[T8l shows the contour plot, for which the azimuth range is di¬ 
vided into bins of 10° and the logarithm of the galactocentric 
radius into bins of 0.2 and the clump mass is summed in each 
bin. The peaks in the mass correlate well with the Scutum-Crux 
arm, Sagittarius and Perseus arms. The map reveals that the mass 
in the spiral arms is clumpy. The distribution of the ATLAS¬ 
GAL sources in the first quadrant is consistent with the 13 CO 
surface brig htness of GRS clouds in th e first quadrant shown 
in Fig. 12 of iRoman-Duval et al.l (2009). In addition, our com¬ 
plexes with the largest number of sources (see Table IB. it . which 
are known molecular cloud complexes, are presented as red dots 
in Fig. [18] and mostly coincide with the enhancements of the 
summed mass distribution. Two complexes, NGC6334 and N49, 
do not correlate well with the mass peaks. They have the smallest 
summed clump mass of the complexes with the largest numbers 
of sources in the first and fourth quadrant, given in Table IB. ll In 
addition, two mass peaks are not associated with any of the com¬ 
plexes containing the largest numbers of clumps from Table lB.il 
One of them is at 6 « -127° and R ga i ~ 4.7 kpc and coincides 
with complex 78 (/ = 18.197°, b = -0.308°), associated with the 
HII region G18.2-0.3. Although it harbours fewer sources, 11 
clumps with observed velocities, it is located at the far distance 
of 11.9 kpc and has a large summed clump mass of 4.64 x 10 5 
M 0 . The other mass enhancement at 0 ~ 110° and R ga i ~ 4.6 
kpc is associated with complex 242 (/ = 336.97°, b = -0.01°), 
which contains 15 sources with observed vel ocities a nd coin¬ 
cides with the HII region G336.84+0.047 dCaswell & Havnesl 
119871) . Our assignement to the far ki nematic distance of 10.9 kpc 
is consistent with p revious studies: iJones & Dickevl (1 201 2l) and 
iGarcfa et all (1201411 also place the HII region at the far distance 
using HI absorption. The far distance leads to a large summed 
clump mass of 5.07 x 10 5 M Q of complex 242. 

The correlation between the location of a large number of 
high mass star forming clumps and spiral arms as provided by a 
four-armed model of the Galaxy (see Fig. IT8l) hints at star for¬ 
mation activity ongoing inside the arms. Because we expect that 
this trend is also reflected by the distribution of the gas mass, we 
investigate if the largest amount of gas mass is confined to spiral 
arms. We therefore study the variation of the gas mass with the 
galactocentric radius in Fig. [19] The upper panel shows the mass 
distribution of the ATLASGAL sources in the first quadrant and 
the lower panel that of the sample in the fourth quadrant plotted 
against galactocentric radii. The black points indicate masses of 
the whole sample above the completeness limit of 1000 M e (see 
Sect. 16.31) and the red points show the summed mass surface den¬ 
sity in each bin of 1 kpc in galactocentric radii with the error (ler) 
calculated from Gaussian error propagation. 

In the first quadrant there is a large number of masses be¬ 
tween ~ 1000 and 10 4 M 0 in addition to large mass values up to 

4 x 10 4 M q at R( ja | ~ 4 kpc resulting in a rising sum of the mass 
surface density per bin to 9 x 10 4 M e kpc -2 , which agrees with 
the increase in the number of ATLASGAL sources at approxi¬ 
mately 4 kpc in Fig. [6] After a decrease to 2.1 x 10 4 M e kpc -2 at 

5 kpc, the summed mass surface density is increasing to 2.5 x 10 4 
M 0 kpc -2 , which is consistent with the distribution of galacto¬ 
centric radii showing another peak at approximately 6 kpc in the 
upper plot of Fig. [6] 

The sum of the mass surface density in the fourth quadrant 
(see lower panel of Fig. IT9l ) is rising to 2.1 x 10 4 M Q kpc -2 
at a galactocentric radius of ~ 3 kpc, which is reflected by 


the slightly increasing number surface density in the lower his¬ 
togram of Fig. [6] After the summed mass surface density is 
slightly decreasing up to 4 kpc, it enhances to 4.6x 10 4 M e kpc -2 
at 5 kpc. This trend is also seen by the galactocentric distribution 
in Sect. 16. II There is a decrease of the summed mass surface den¬ 
sity at galactocentric radii larger than 6 kpc with a slight rise at 
7.5 kpc, which is consistent with the distribution of the number 
surface density. This analysis consequently reveals an agreement 
of the gas mass with spiral arms, which are also correlated with 
the number of ATLASGAL clumps as shown in Sect. 16.11 

We also investigate the variation of the mean gas mass es¬ 
timated in bins of 1 kpc in galactocentric radii. However, the 
peaks of the surface density are not revealed by the average 
gas mass, which is constant within the errors within the in¬ 
ner Galaxy. This reveals that the clumps in the spiral arms are 
not more massive, but that the enhancement of the mass in the 
spir al arms results fro m a larger number of clum ps. In addi¬ 
tion, lEden et all d2012l) and lBattisti & Heveil (12014l) revealed no 
changes in the clump formation efficiency with different envi¬ 
ronments within the Galactic plane. A recent study of differences 
in the clump formation efficiency along the line of sight within 
37.83° < / < 42.5° (lEden et al.112013b results in no variation 
in the clump formation efficiency between arm and interam re¬ 
gions. The rise of the mass in the spiral arms is therefore not 
caused by an increased clump formation efficiency, but likely by 
a flow of the gas of the interstellar medium and molecular clouds 
in high density regions resulting in source crowding and an in¬ 
creased probability of collisions in the spiral arms. 

7.4. Clump mass function 

7.4.1. Mass distribution of ATLASGAL clumps and clusters 

The derivation of the power-law exponent of the clump mass 
function (see Sect. I67fl allows us to analyse if it is consistent 
with the stellar initial mass function (IMF). This would favour 
theoretical models, which explain high mass star formation as 
extension of the theory of low mass stars forming from a molec¬ 
ular cloud collapse via disk accretion. Our a average of -2.29 
for the differential mass functions agrees w ith the power-law ex¬ 
ponent of the stellar IMF given bv lKrou pal (120081) of -2.3 ± 0.5. 
This hints at an origin of the IMF from the clump mass function. 
However, the confirmation of models of high mass star formation 
requires that clumps are converted one-to-one into stars or bina¬ 
ries, which cannot be revealed from our data because one clump 
is likely to form a whole cluster of stars. The similarity between 
the slope of the clump mass function and the IMF shows that the 
star formation efficiency in clumps does not depend on the mass 
of the clump. 

Studies of embedded young clusters iLada & Ladai I2003I) . 
young and old globular clusters, massive open clu sters , and as- 
socia t ions as well as interstellar c louds (lElmegreen & Efremov! 
1 19971 : Portegies Zwart et al.l2010ll yield a power-law exponent a 
of ~ -2, which lies in the a range determined from our clump 
mass function fits. This consequently reveals a consistency of 
the mass function of clusters and ATLASGAL clumps. The mass 
spectrum is therefore not affected by fragmentation. This result 
is also obtained by other studies using d ust c ont inuum emis- 
sion as well as of high-density gas tracers. iBeltran et ail ( 2006 ) 
measured the 1.2 mm dust continuum of clumps in high mass 
star forming regions; a fit to the mass spectrum above 92 M e 
of sources within 6 kpc yields a slope of -2.1 + 0.02, which is 
similar to that derived for pre- and protocluster masses and also 
agrees with our a values of the ATLASGAL subsample from 2 to 
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Fig. 18: The logarithm of the galactocentric radius plotted against the azimuth, 6, of the ATLASGAL sources in the first and fourth 
quadrant above 1000 M e as scatter plot in the left panel, clumps within 2 kpc of the galactic centre are excluded. The angle 6 to 
the Galactic centre-Sun axis is counted anticlockwise as viewed from the north Galactic pole. The right panel displays the contour 
plot, for which we summed the clump mass in each azimuth bin of 10° and each logarithmic galactocentric radius bin of 0.2. The 
contours give 10 to 90% in steps of 10% of the peak summed clump mass per bin. The spiral arms originate from the model by 
IValleel (Il995l) . The black, red, green, and blue lines indicate the Sagittarius arm, the Perseus arm, the 3-kpc arm, and the Scutum- 
Crux arm. Complexes with the largest number of sources are illustrated as red dots, most of them agree well with the enhancements 
of the summed mass. 


5 kpc. lShirlev et al. ( 2003 ) observed high mass cores, which will 
form clusters or associations, in the high-density probe CS and 
calculated virial masses. Fitting the cumulative mass function to 
their distribution above 1000 M 0 gives a power-law exponent of 
-1.91 + 0.17, slightly smaller than our a value while still simi¬ 
lar to our power-law exponent and also co nsistent with the IMF 
slope of stars within clusters (tShirlev et al.112003 0. This again, 
shows the similarity of the clump and cluster mass function. 


7.4.2. Comparison with literature values 


iReid & Wilsonl ( 20061) investigated the mass distribution of star 
forming regions with different median masses. They fitted dou¬ 
ble power laws to the cumulative mass functions of clumps or 
cores in low mass objects such as found in the p Oph region, 
as well as in complexes of high mass star formation, e.g. W43 
and RCW 106. Their (Thigh values lie between -1.8 and -3.1 
with a mean of -2.4 + 0.1. They find no dependence of (Thigh on 
the median clump mass for a mass range from 0.4 to 1000 M 0 . 
The ATLASGAL sample covers the same mass range as well 
as sources of still higher masses, from 0.4 M 0 up to 5.6 x 10 4 
M 0 . Fits to its differential and cumulative mass functions yield 
power-law exponents, which are consistent with the mean given 
by iReid & Wilson ( 20061) . 

We also study the dependence of the power-law exponents 
on the mass obtained by more recent studies of high mass star 
forming regions. Table 0 gives the name of the source sample, 
the mass range, the distance range, the minimum mass, A7 mm , 
above which the distribution is fitted, the power-law exponent, 
and the reference. We determine the minimum mass from the 
differential mass functions shown in the listed articles. We con¬ 
sider source samp le s cove ring different mass ranges from cores 
in IRDCs dRathborn e et al .11 20061) wit h masses up to ~ 2000 M 0 , 
clumps in IRDCs (iRagan et al.l 20091) exhibiting up to 3000 M 0 , 
high mass clumps detected in the 1.2 mm dust continuum with 
masses up to 4100 M 0 (Beltran et al.l2006l) . clumps in molecular 


cloud complexes from the GRS (Simon et al .11200 1 1) having up to 
10 4 M 0 as well as IRDCs (Peretto & Fulledfcoiol) with masses 
up to 10 5 M 0 . Comparing the different mass distributions reveals 
that the IRDCs and the ATLASGAL sources have the broadest 
mass ranges with a similar minimum mass. We plot the power- 
law exponents of the objects listed in Table [5] against their min¬ 
imum masses in Fig. [20] The study of GRS molecular clouds 
(see Table [3 with varying star formation activity and located at 
various distances dSimon et al.l 1200 II) revealed no difference in 
the slope of the mass functions. Figure l20lshows that it is also in 
agreement with the power-law exponent of clumps in IRDCs and 
the ATLASGAL sample, while the slope of the mass spectrum of 
1.2 mm dust clumps and of cores in IRDCs are slightly steeper. 
This comparison indicates that the mass distribution of cores, 
clumps, and clusters consisting of several clumps can be fitted 
by power laws with similar slo pes. We also add four high mass 
star forming regions studied bv lSwift & Beaumont! ( 2010 ). who 
estimated the slope of th e mass spectra u sing a method based on 
maximum likelihood (IClauset et al.l12009!) . These samples con¬ 
sist of M17, NGC 7538, Cygnus X, and NGC 6334, which are 
all located nearby and exhibit various mass ranges: Sources in 
M17 only have masses up to 120 M 0 , while Cygnus X con¬ 
tains clumps with up to 949 M 0 , still larger mass ranges are 
probed by NGC 7538 exhibiting masses up to 2700 M 0 an d NG C 
6334 with masses up to 6000 M 0 . Although ISwift & Beaumontl 
(2010 ) used a method, which is independent of the histogram 
binning, their power-law exponents are similar to those of the 
other samples in Table [5] The slope of NGC 6334 is consistent 
with that of clumps in IRDCs, of ATLASGAL sources, and of 
clumps in GRS molecular clouds. The power-law index of M17, 
NGC 7538, and Cygnus X are the same as obtained for cores in 
IRDCs and of 1.2 mm dust clumps. Figure l20l shows that there 
is no trend of a with the minimum mass. Table [^reveals that the 
power-law exponents of the mass function derived from 13 CO 
emission, observed towards GRS molecular clouds, and from 
dust emission are consistent. 
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Fig. 19: Distribution of gas masses above 1000 M Q is illustrated 
as black points in the upper panel for the ATLASGAL sources in 
the first quadrant and in the lower panel for clumps in the fourth 
quadrant. Red points indicate the sum of the gas mass surface 
density estimated in bins of 1 kpc in galactocentric radii. 


Using CO emis sion of molecular clouds in the inner Galaxy 
iRosolowskvl (120051) fitted a slope of -1.5 to the mass spectrum, 
which is flatter than our range of a values. This difference re¬ 
sults from a larger amount of material probed by CO, which 
traces low-density regions. High-density tracers such as NH 3 de¬ 
tected in the ATLASGAL sample are more appropriate to study 
the clump mass distribution because stars form only in molecular 
clouds that exceed a gas surface density threshold of I gas ~ 116 
M g pc -2 (lLada et al.]l2012l) . 

To investigate any differences between the mass distributions 
of sources in different evolutionary phases, we distinguish be¬ 
tween clumps that have the ratio of the MSX 21 pm flux to the 
870 pm peak flux density smaller or larger than 0.8 as described 
in Sect. 16.51 This leads to a cold sample with a mean gas mass 
of 2980 M g and a warm sample with a mean mass of 3325 M 0 . 
Figure [21] compares the mass spectrum of the cold sample in 
red with that of the whole ATLASGAL sample in black, which 
shows that the two have the same distribution. From a power-law 
fit to the differential clump mass function of the whole ATLAS¬ 
GAL sample above the completeness limit of 1000 M 0 we derive 
a slope of -1.83 ± 0.1. The fit to the mass spectrum of the whole 
ATLASGAL sample in Fig. |2T1 indicates that it has a lognormal 
distribution (see Sect. PI We obtain a power-law exponent of 
-1.86 + 0.09 from the fit of the mass distribution of the cold 


sources above 1000 M 0 , which agrees within the errors with the 
slope of the whole ATLASGAL sample. The clump mass distri¬ 
bution is therefore already set in IRDCs and does not change in 
lat er evolutionary phases . This is also supported by the analysis 
of lUrquhart et al I (12013al) . who derived a similar power-law ex¬ 
ponent of -2.0 ±0.1 for ATLASGAL sources associated with 
methanol masers. A comparison of their clumps with other sam¬ 
ples also yields a similar clump mass function for the embedded 
phases of high mass star formation. lOlmi et akl (12013l) divided 
Hi-GAL sources at a longitude of 30° and 59° in samples of star¬ 
less and proto-stellar clumps. They fitted the differential clump 
mass function to the distribution of the different samples and also 
obtained a similar power-law index for the whole sample, the 
starless clumps, and proto-stellar clumps. In addition, we com¬ 
pare our result with the mass spectra of mas sive dust clumps in 
the fourth quadrant from lBeltran et al.l t 2 006 ) who investigated if 
these are associated with point-like, diffuse or no MSX emission. 
Their study of different clump properties revealed that their sub¬ 
sample without MSX counterpart has a lower mean mass of 96 
M 0 than sources associated with point-like or diffuse MSX emis¬ 
sion, which has a mean mass of 336 M 0 similar to the trend of 
average masses of our subsamples. We cannot fit the mass distri¬ 
bution of clumps not as socia ted with mid-IR emission from the 
sample of lBeltran et al.l (120061) because of the small source num¬ 
ber, which shows that large samples such as ATLASGAL are 
required for a statistically significant analysis. The power-law 
fit to the differential mass distr ibuti on of massiv e dust clumps 
with MSX counterparts from lBeltran et al.l (2006j) gives a slope 
of -1.84 ± 0.06, which is consistent within the errors with our fit 
result of the whole ATLASGAL sample. 


7.4.3. Variation of the clump mass function with 
Galactocentric radii 

To investigate the change in the slope of the mass distribution 
in the inner Galaxy we divide ATLASGAL sources with galac¬ 
tocentric radii between 3 and 8 kpc into equally spaced bins of 
1 kpc and derive their mass spectra above 1000 M 0 , which is 
the turnover in the complete mass distribution as seen in Fig. [9] 
Fitting of differential clump mass functions results in power-law 
exponents between -1.32 and -1.86. They are plotted against 
galactocentric radii in Fig. [22] which shows a slightly decreasing 
slope with increasing galactocentric radius out to ~ 7 kpc after 
which the slope increases to a = -1.3 at 8 kpc. This indicates 
that there is generally a higher proportion of lower mass clumps 
as a function of galactocentric radius. To compare the mass 
function of inner and outer Galaxy data we use o bservatio ns 
of 132 massive protoclusters in the outer Galaxy (iKlein et al.l 
l2005[) . Their gas masses range from ~ 0.2 to 5200 M 0 with an 
average of 290 M 0 . We fit a power law to the differential mass 
spectrum of the protoclusters and derive an exponent of -1.48, 
which is flatter than measured for the ATLASGAL sample. 
However, the outer Galaxy sample consists of fewer sources 
than the observed ATLASGAL clumps, which results in a larger 
uncertainty of the protocluster mass spectrum. 


7.5. How many ATLASGAL sources form high mass stars? 

The derivation of the kinematic distances allows us to examine 
the range of gas masses and sizes of the ATLASGAL sample 
(see Sect. IQ) . which consists of different kinds of sources. To 
examine their star formation activity we have to connect masses 
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Table 5: Fit results of the mass spectra of different high mass star forming samples. 


Sample 

mass range (M 0 ) 

distance range (kpc) 

Mmin 1 (Ms) 

a 

Reference 

IRDC cores 

10-2100 

0.4 - 7.8 

100 

-2.1 + 0.4 

Rathborne et al. (2006) 

IRDC clumps 

30 - 3000 

2-5 

40 2 

-1.76 + 0.05 

Ragan et al. (2009) 

IRDCs 

O 

1 

o 

Ul 

2.5 - 7.5 

800 

-1.85 + 0.07 

Peretto & Fuller (2010) 

1.2 mm dust clumps 

10-4100 

0.01-4 

92 

-2.1+0.02 

Beltran et al. (2006) 

GRSMC 45.60+0.30 clumps 

4 - 1000 

1.8 

15 

-1.78 + 0.13 

Simon et al. (2001) 

GRSMC 43.30-0.33 clumps 

10-1600 

3 

40 

-1.85 + 0.12 

Simon et al. (200 ^ 

GRSMC 45.46+0.05 clumps 

25 - 10 4 

6 

250 

-1.80 + 0.10 

Simon et al. (200D 

W49 

250 - 4 x 10 4 

11.4 

1000 

-1.52 + 0.20 

Simon et al. (200") 

ATLASGAL 

0.4 - 5.6 x 10 4 

0.2 - 17 

1000 

-1.83 + 0.1 

Sect. |6.4| 

M17 

0.8 - 120 

1.6 

3 3 

-2.1 + 0.2 

Swift & Beaumont (2010) 

NGC 7538 

1.4-2700 

2.8 

15 3 

-2.1+0.36 

Swift & Beaumont (2010) 

Cygnus X 

4-949 

1.7 

5.3 3 

-2.1 + 0.2 

Swift & Beaumont (2010) 

NGC 6334 

3 - 6000 

1.7 

32 3 

-1.7 + 0.1 

Swift & Beaumont (2010) 


Notes. (i) M m j n is the minimum mass, abov e which the distribution is fitted. 

We u se the mass mum giv en in Table 1 in lRasan et al.! 1200 9.1 as the minimum mass. 

(j hSwift & Beaumont! d2010h fit a double power law to the mass spectrum of IRDC clumps, we use the break mass as the minimum mass. 



Fig. 20: Slope of the mass functions is compared with the mini¬ 
mum mass for different high mass star forming samples from Ta¬ 
ble 0 Cores within IRDCs are indicated as black point, clumps 
within IRDCs as purple point, IRDCs as green point, 1.2 mm 
dust clumps as blue point, clumps within GRS molecular clouds 
and W49 as light blue points, ATLASGAL sources as black tri¬ 
angle, M17, NGC 7538, Cygnus X, and NGC 6334 as red trian¬ 
gles. 


and sizes, rela tions between the two h ave been obtained by pre¬ 
vious studies. iHeiderman et al.l (1201 Ot) fitted the star formation 
rate against surface densities of 20 molecular clouds, associated 
mainly with low mass star formation, as well as high mass dense 
clumps. The changing slope of the fit from a steep to a linear 
function yields a star formation threshold of 129 ± 14 M 0 pc 2 . 
lLada et all ( 2010 ) analysed the relation between the star forma¬ 
tion rate of molecular clouds within 500 pc and the mass of the 
cloud. This revealed that the star formation activity does not de¬ 
pend on the total mass of the cloud, but on the extinction of the 
molecular gas and therefore volume density. lLada et aTI (1201 Oi) 
obtained a linear relation between the star formation rate and the 
mass of a cloud above a threshold in visual extinction of ~ 7 
mag, which is a gas surface density threshold of ~ 116 M e pc~ 2 . 
The analysis of molecular clouds, which do not form high mass 


Fig. 21: Similar slopes result from power-law fits to the differen¬ 
tial mass function of the whole ATLASGAL sample in black and 
of a cold subsample in red with the ratio of the MSX 21 fjm flux 
to the 870 fim peak flux density smaller than 0.8. The vertical 
line displays the completeness limit of 1000 M 0 , above which 
the mass spectra are fitted. 


stars such as Taurus, Pe rseus, Ophiuchus, the Pipe Nebula, by 
[Kauffmann et al.1 (1201 Oi) revealed that their masses are limited 
by 

M(r) = SSOMoCReffpc- 1 ) 1 ' 33 (9) 


with the gas mass, M(r) 
given in 


and the eff ectiv e radius, /Air. as 


uiv v ii i-j. r w j. ^ 

Contreras et al.l (1 201 3l) and lUrauhart et al.l (1201 4al) . 
Kauffmann et al.1 (1201 01) de creased the dust opacities of 


Ossenkonf & Henning! (1 1 9941) by a factor of 1.5, while we use 


the unmodified opacities in the calculation of the gas mass. 
This change is taken into account in equatio n |9l where we re¬ 
duce the factor 870 given in lKauffmann et all (1201 Ot) to 580. In 
contrast to the nearby clouds, two massive star forming com- 
plexes have masses larg er than derived in equation [9] which led 
iKauffmann et al.l (1201 0i) to consider the relation as threshold for 
high mass star formation. However, this must still be confirmed 
by the study of larger samples of clouds that form high mass 
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Fig. 22: Fitting of the differential mass distributions of ATLAS- 
GAL clumps above 1000 M Q within 1 kpc over the range of 
galactocentric radii from 3 to 8 kpc gives similar power-law ex¬ 
ponents, a. There is a slightly decreasing slope of a with increas¬ 
ing galactocentric radius out to ~ 7 kpc and an increasing slope 
to a = -1.3 at 8 kpc. 

stars. R ecent work bv lUrauhart et alJ (l2013al) and lUrauhart et alJ 
(1201 3bl) has analysed the mass-size relation of a large number of 
ATLASGAL sources, but it traces only the more evolved phases 
of massive star formation from methanol masers to UCHIIRs. 
In contrast, ATLASGAL clumps studied in this article are se¬ 
lected from the unbiased ATLASGAL sample and therefore in¬ 
clude cold sources in an early, even prestellar, evolutionary stage 
as well as objects in later phases. 

We show the correlation between the logarithm of the ef¬ 
fective radius and of the gas mass of the ATLASGAL sam¬ 
ple with derived kinematic distances in Fig. [23] The light blue 
line shows the l imits in gas surface density of 116 M 0 pc 2 
(iLada et al.l2010t) and 129+ 14M 0 pc" 2 (iHeiderman et al.l2010l) . 
above which sources are likely to form stars. Using an average 
of the two thresholds, 122.5 M 0 pc ~ 2 , we obtain 3226 ATLAS¬ 
GAL clumps with masses larger than this limit, which is 100% 
of the sources with determined gas masses. The comparison to 
BGPS sources dDunham et al.ll20l71) reveals that a smaller num- 
ber of them, ~ 46 % of the whole sample , satisf y the criteria from 
ILada et al.l d20loh and IHeiderman et akl d2010l) using the surface 
density of the clumps. However, this fraction will increase to 
~ 87% if they refer to smaller regions within the sources by tak¬ 
ing the peak column density into account. 

The red line i n Fig. 1231 indicates the mass-radius relation¬ 
ship derived by iKauffmann et al.l d2010h . We obtain 2953 AT¬ 
LASGAL sources, 92% of the whole sample, which exhibit 
larger masses than expected from equation[9]and therefore might 
form high mass protostars. This fraction of ATLASGAL clumps 
is higher than the result from iDunham et al.l d201 ll) . who cal¬ 
culated that ~ 48 % of their sample have masses larger than 
the limit given by IKauffmann et al.l ( 2010 ). and smaller than 
~ 97% of ATLASGA L sources identified as methanol masers 
dUrauhart et al.112013a ) that are potentially forming high mass 
protostars. While those trace later evolutionary phases from hot 
molecular cores to UCHIIRs, the ATLASGAL clumps anal¬ 
ysed in this article and the BGPS sources are unbiased sam¬ 
ples covering various evolutionary stages. Because these also 
include cold sources in an early phase, which show no sign of 
high mass star formation yet, they can explain the larger por¬ 


tion of clumps below the threshold from iKauffmann et al.l 1 2010 1 
compar ed to ATLASGAL so urces associated with methanol 
masers. lUrauhart et al.l d2013al) obtain similar results from a fit 
to the mass-size relation of ATLASGAL sources associated with 
methanol masers and with compact HII regions. A fit to our data 
yields logio(M) = 3.22 ± 0.002 + (1.76 ± 0.01) x logio(R e ff) i n 
Fig. [23] Their comparison shows that the slope of our sample, 
which also contains cold clumps, is consistent within 3cr with 
that of the sources in later evolutionary phases, although Fig. [23] 
reveals the trend of a broader mass range of the ATLASGAL 
sample studied in this article extending down to lower values. 
The dashed black line in Fig. [23] which indicates the sensitiv¬ 
ity of the ATLASGAL survey of 5cr, is a lower boundary to the 
derived gas masses ._ 

ICsengeri et al.l (120141) used a limit of 650 M e , above which 
ATLASGAL sources are likely to harbour massive dense cores 
and high-mass protostars. This is indicated as a dotted green 
line in Fig. We find 2097 ATLASGAL clumps, 71% of 
the sources with mass es larger than the threshold given by 
IKauffmann et al.l d2010t) . that also exhibit masses larger than 650 
M 0 . 

The dotted blue line in Fig. 1231 shows masses corresponding 
to a constant surface density of 1 g cm" 2 . This value is deter¬ 
mined to be the limit abov e which fragmentation is su ppressed 
and high mass stars form dKrumholz & McKcejl2008 ). The in¬ 
vestigation of high mass star forming regions as well as star clus¬ 
ters with massive stars bv lMcKee & Tanl (120031) revealed that the 
surface densities of ~ 1 g cm " 2 of these objects vary by a factor 
of 4. Because this threshold is derived for small-scale structures 
such as cores, while we are tracing larger objects with clump 
radii and smaller densities, only a small part of the ATLASGAL 
sources lie close to the dotted blue line. Details about clumps, 
which have high gas masses concentrated to small source radii, 
are given in Appendix lB.il 

To analyse if there is a trend of star formation activity with 
the mass-size relation of ATLASGAL clumps with derived kine¬ 
matic distances we divide them into two samples with the ratio 
of the MSX 21 pm flux to the 870 pm peak flux density smaller 
or larger than 0.8. As described in Sect. 17.4| this provides us a 
cold and a warm sample. We plot the logarithm of their gas mass 
against the logarithm of the effective radius as red triangles and 
black points in Fig. [24] The two subsamples exhibit similar dis¬ 
tributions. The lines illustrate the same relations as given in Fig. 
1231 In addition, we average the gas mass within bins of 0.5 in the 
logarithm of the effective radii to examine if there is a trend in 
the mean values, which are shown in the lower panel of Fig. [24] 
However, we find no difference between the distributions of the 
two subsamples either. 

In conclusion, most ATLASGAL clumps exhibit gas masses 
larger than the limit above which high mass stars form 
(iKauffmann et al.!l201Q ). This indicates that the majority of AT¬ 
LASGAL clumps have the potential to form high mass proto¬ 
stars. In addition, the distribution of the mass-size relation is sim¬ 
ilar for a cold subsample in an early evolutionary stage without 
star formation activity as well as for warm clumps harbouring an 
embedded heating source. 

8. Conclusions 

We calculated near and far kinematic distances from the radial 
velocity of various molecular lines, mostly the NH 3 ( 1 , 1 ) inver¬ 
sion transition, to a large sample of dust c lumps discovered by 
the ATLASGAL survey using the Brand & Blitzl (1 19931) rotation 
curve. This section summarizes our results. 


Article number, page 24 of[36] 



































































M. Wienen et al.: ATLASGAL - Kinematic distances and the dense mass distribution of the inner Galaxy 


o 

CT> 

O 



- 2-10 1 - 2-10 1 


og 1 0(Effective Radius/pc) 


og 10(Effective Radius/pc) 




log 1 0(Effective Radius/pc) 


Fig. 23: Logarithm of the gas mass compared with the log¬ 
arithm of the effective radius of ATLASGAL sources. The 
red line illu strate s the threshold of massive star formation by 
iKauffmann et al.l (l2010h . The dotted green line illustrates a 
threshold of 650 M Q for the formation of high mass dense cores 
and massive protostars derived by ICsengeri et al.l (120141). Limit 
in ga s surface density as a verage of 116 M r . ; pc 2 (lLada et al.l 
[ 20id ) and 129+ 14M 0 pc~ 2 dHeiderman et al.lBoiOt) is displayed 
as light blue line. A gas surface density of 1 g cm -2 is shown as 
dotted blue line. The sensitivity of the ATLASGAL survey is 
displayed as dashed black line. The upper panel illustrates the 
contour plot with a binning of the logarithmic effective radius of 
0.2 and the logarithmic gas mass of 0.2. 


1. To determine distances we first divide the ATLASGAL sam¬ 
ple into groups of sources, which are coherent in space and 
velocity. We associate observed ATLASGAL sources with 
molecular clouds by comparing their location with 12 CO 
emission probing molecular clouds on a large scale, with the 
13 CO (1-0) line in the first quadrant and 870 fim dust con¬ 
tinuum tracing the small-scale structure. This allows us to 
assign additional ATLASGAL sources without any known 
velocities to the complexes of measured clumps and there¬ 
fore the same distance. 

2. Near and far kinematic distances to the complexes of AT¬ 
LASGAL _sourcesjtre computed using the rotation curve by 
iBrand & Blitzi d 1993b . To resolve the KDA we investigate HI 
self-absorption toward all ATLASGAL clumps and HI ab¬ 
sorption toward 65 sources in the first quadrant with an em- 



Fig. 24: The mass-size relation for a cold and a warm ATLAS¬ 
GAL subsample (see Sect. 17.51 ) is shown as a red and black con¬ 
tour plot in the top panel, bins of 0.25 are used for the logarithm 
of the effective radius and bins of 0.2 for the logarithm of the gas 
mass. The scatter plot with the cold subsample as red triangles 
and the warm sources as black points is displayed in the middle 
panel. Mean values of the logarithm of the gas masses and of the 
effective radii are plotted in the bottom panel. The lines are the 
same as described in Sect. 17.51 


bedded HII region and 146 sources in the fourth quadrant. 
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We combine our KDA resolutions derived from the two tech¬ 
niques and obtain distances to 689 complexes. We locate 143 
complexes in the first quadrant at the near distance, 148 at the 
far distance, and 5 at the tangent point, while 234 complexes 
in the fourth quadrant are assigned to the near distance, 141 
to the far distance, and 18 at the tangent point. 

3. While Galactocentric distances of the ATLASGAL subsam¬ 
ple in the first quadrant show a concentration of sources at 
~ 4 and 6 kpc, Galactocentric distances in the fourth quad¬ 
rant reveal an increase of sources at 3, between 5 and 6 kpc, 
and ~ 7 kpc. These enhancements of clumps probe different 
spiral arms of the Milky Way. 

4. An exponential fit to the number distribution of ATLAS¬ 
GAL sources with the distance from the Galactic plane 
gives a similar scale height of ~ 28 ± 2 pc for the sam¬ 
ples in the first and fourth quadrant. In addition, we ob¬ 
tain a similar displacement of ATLASGAL sources below 
the Galactic midplane in the first and fourth quadrants, at 
~ —10 ± 0.5 and ~ -4+1.7 pc. This is consistent with 
the r esults derived for other high mass star forming sam - 
ples dGreen & McClure-Griffithsl[201 It iBecker et al .11 1 99(3) . 
The distribution of the scale height and the distance from the 
midplane are constant from 3 kpc to 8 kpc with respect to the 
outer Galaxy. 

5. Using kinematic distances we calculate a broad mass and 
size range of the ATLASGAL sources. Within 2 to 18 kpc we 
trace clumps with 50 to 500 M Q and radii from 0.15 to 1.5 pc. 
A fit of the differential clump mass function to the mass dis¬ 
tribution of ATLASGAL sources between 2 and 5 kpc results 
in a power-law exponent of-2.2+0.1 on average. The cumu¬ 
lative mass function of ATLASGAL clumps is better repre¬ 
sented by a lognormal distribution than by a double power 
law. Our slope of the differential mass dis tributi on agree s 
with the power-law exponent of the IMF (iKrounal 120081) . 
the mass function of clusters, associations, and that derived 
from high-density tracers as well as d ust c ontin uum emis¬ 
sion for massive star forming regions (lLada & Ladal l2003t 
iBeltran et alJ 120061 : ISimon et alj|200ll) . Our fit to the mass 
spectrum of different ATLASGAL subsamples reveals that 
the power-law index does not vary from the early, cold phase 
to later evolutionary stages. In addition, the slope also re¬ 
mains similar within the inner Galaxy. 

6. The comparison of mass and radius shows that the 
whole ATLASGAL s amp le h as gas surface densities higher 


_ ig< _, 

than 116 M f .,pc 2 ( Lada et all 1201 Oh and 129 M 0 pc 2 
dHeiderman etal.li201 (?) . which is the threshold for star for¬ 
mation. Masses of 92% of the ATLASGAL clumps are 
larger than the limit for high mass star formation given by 
iKauffmann et al.l (l2010t) . Most sources might therefore po¬ 
tentially form high mass protostars. A fit of the mass-size 
relationship by a power law results in a slope of 1.76 ± 0.01, 
which remains invariant for subsamples in different evolu¬ 
tionary phases. 
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Appendix A: Grouping of ATLASGAL sources 

We give in this section details about our analysis to group 
sources, which are coherent in space and velocity us¬ 


ing the friends-of-friends algorithm dHuchra & Gelled 119821 : 
iMoore et al.lll993l : iBerlind et al.ll2006h . As a first step we sort 
ATLASGAL sources with known velocities (see Sect. IT2l> in 
the first and fourth quadrant according to an increasing galac¬ 
tic longitude range. For each ATLASGAL source we search for 
related sources in the whole sample according to a predeter¬ 
mined threshold in position and velocity. Associated ATLAS¬ 
GAL sources are added to the same group until no more sources 
that lie within a maximum distance and velocity interval to the 
neighbouring sources are found. The group is then considered 
to be complete and a new search around the next ATLASGAL 
source of the catalogue begins. All distances between at least two 
individual sources within one group must not exceed the maxi¬ 
mum distance and the velocities of these must be within a given 
velocity interval. However, the size and the maximum difference 
of velocities within a group can exceed the input values for the 
maximum distance and the velocity interval used to search for 
associated sources. This algorithm results in an arrangement of 
the ATLASGAL sample into groups and sources that do not be¬ 
long to any group. 

We study the dependency of the grouping on the input pa¬ 
rameters by varying the maximum distance and the velocity 
interval. This results in various divisions of the ATLASGAL 
sources into groups, which are analysed to improve the match¬ 
ing between our grouping and properties of known molecular 
cloud complexes (see Appendix IB. 2 k For different input values 
for the maximum distance, A d u and a fixed interval of 10 km s 1 
for velocities of neighbouring sources in a group we plot his¬ 
tograms of the size of a group in Fig. lA. II The size is defined as 
(A/ x Ab ) 1 /2 with the maximum difference in longitude. A/, and 
in latitude, Ab, of sources within a group. For Ac/; = 0.05° or 
0.1° the size of most groups is smaller than those values, while 
the size of a larger number of groups is in the range of the input 
value Adi = 0.2° or 0.3°. If Adi is larger, 0.5°, some groups will 
have sizes larger than Ac/;. This trend is visible for the ATLAS¬ 
GAL samples in the first (upper panel of Fig. IA.1I) and fourth 
quadrant (lower panel of Fig. IA. II) . although it is more obvious 
for the groups in the fourth quadrant. To avoid gathering sep¬ 
arate groups together, we do not take a large input value such 
as 0.5°, but choos e 0,3°. In addition, the m edian size of GRS 
molecular clouds (iRoman-Duval et all l2009t) of ~ 0.3° is also 
consistent with the maximum distance of ATLASGAL sources 
within a group, which corresponds to ~ 22 pc at a m edian dis¬ 
tance of ~ 4 kpc. Using 12 CO emission iGarcfa et al.1 (120141) ob¬ 
tained even sizes up to 100 pc for giant molecular clouds in the 
fourth quadrant. 

Figure [701 illustrates histograms of the maximum difference 
in velocities within a group for varying input values of the ve¬ 
locity interval used to search for neighbouring sources within a 
group, Av;, and a fixed maximum distances between individual 
sources in a group. Figure IA.2I indicates that with an increase 
in Av; from 10 to 15 km s _1 the maximum velocity difference 
in a few groups becomes greater than the chosen Avj. This is 
similar to the trend seen in the maximum distance, but less ob¬ 
vious. Since the maximum velocity difference of most groups 
should not exceed the input value of Av;, we take 10 km s' 1 for 
that. This analysis lets us choose the arrangement of ATLAS¬ 
GAL sources into groups resulting from 0.3° and 10 km s _1 as 
input values for the maximum distance and the maximum dif¬ 
ference in velocities of neighbouring sources in a group. The 
maximum velocity difference of ATLASGAL sources in a group 
also agrees with the velocit y dispersion of GRS molecular clouds 
dRoman-Duval et al.l2009l) ranging up to 10 km s -1 and with the 
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median linewidth of GMCs in the fourth quadrant dGarcfa et alJ 
120141) of 9 km s~ 1 . 

We call sources that are associated by their position and ve¬ 
locity ”a group” and refer to the group’s spatial structure in 
a comparison with the ATLASGAL dust continuum emission 
as ”a complex”. We illustrate an example of a complex within 
22.1° < l < 22.6° and 0.25° < b < 0.5°, which results from our 
method, in the top panel of Fig. [2] 

This method results in 296 groups in the first quadrant and 
393 groups in the fourth quadrant, an uncertainty associated 
with these identifications might result from gathering separated 
groups along the line of sight together into one group. 

Appendix A.1: Properties of groups 

To characterize the groups we study some physical properties 
such as the velocity dispersion and physical radius. The distri¬ 
bution of the number of sources per group is illustrated in panel 

a) of Fig. IA.3I the histogram of the velocity dispersion in panel 

b) and of the physical radius of the groups in panel c) using dis¬ 
tances determined in Sect. [5] We obtain ~ 300 sources in the first 
and fourth quadrants that within 0.3° and 10 km s' 1 are not re¬ 
lated to another ATLASGAL source, as indicated in panel a) of 
Fig. IA.3I The histograms of the velocity dispersion and physi¬ 
cal radius of the groups containing > 2 sources show that most 
groups exhibit a velocity dispersion of < 2 km s -1 and radius of 
~ 5 pc. The peak in the radius distribution mainly results from a 
subsample of about 100 groups, which contain only two sources 
per group. In contrast, the peak of the velocity dispersion is a 
property of the whole sample. Figure IA.3I re ve a I s a large disper¬ 
sion of radii, which range to 80 pc. We obtain 12 groups out of 
the whole sample of 689 with radii larger than 35 pc, most of 
them consist of more than ten sources. Some of the large radii 
belong to groups at the far distances. 

The correlation of the physical radius and velocity dispersion 
of all 709 groups is shown in Fig. IA.4l A power-law relationship 
between the two prope rties is indicated as black line, which is 
given bv lLarsonl ( ll981 ! 
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6v = 1.43R 0 ' 38 , (A.l) 

where 5v is the velocity dispersion and R the radius. Some of the 
groups are located close to the relationship, but many groups lie 
in a broad distribution. To describe the whole sample we added 
a relation for a gravitationally bound cloud, resulting from a re¬ 
examination of Larson’s law bv lBallesteros-Paredes et ahl (t 2011 ) 

6v = V2 GI.R (A.2) 

with the FL column density I. We plot Eq. |A.2l for FL column 
densities of 10 cm as green line, 10 cm as red line, and 
5 x 10 22 crrT 2 as blue line in Fig. IA.4I The distribution of ra¬ 
dius and velocity dispersion can be fitted by a power law with 
an exponent of 0.4 ± 0.02, which is shown as purple line. The 
slope is thus consistent with the power-law exponent given in 
lLarsoni ( 198 I T We performed a Spearman correlation test for an 
ATLASGAL subsample with heliocentric distances between 2 
and 5 kpc, which gives a correlation coefficient of 0.57 with a 
p-value < 0.0013. This indicates that the correlation between ra¬ 
dius and velocity dispersion is significant over 3cr. A few groups 
are below the green line and exhibit extremely small velocity 
dispersions, <0.1 km s -1 , because they mostly include only 
fewer than five sources with observed velocities. Because these 
do not have good statistics, we aim at describing groups asso¬ 
ciated with a larger number of measured velocities. A contour 


Fig. A. 1: Relative number distribution of ATLASGAL groups 
with the size of a group for different input values of the maxi¬ 
mum distance between individual sources within a group. Ad ;. 
This is varied between 0.05° and 0.5° with a fixed interval for 
velocities of neighbouring sources in a group, Av;, of 10 km s _I 
for the sample in the first quadrant in panel a) and the groups in 
the fourth quadrant in panel b). 

plot is presented in the upper panel of Fig. IA.41 with a binning 
of 0.3 for the logarithm of the radius and for the logarithm of 
the velocity dispersion. It shows a trend of growing radii with 
inc reasi ng vel ocity dispersionsjtnd most groups follow the law 
of lBallesteros-Paredes et all (1201 111 . 

Appendix B: Notes on individual sources 

Appendix B. 1: Sources with extreme properties 
Clumps with high masses 

The comparison of the kinematic distance with the gas mass 
in Fig. [8] reveals some extreme sources clustering at different 
distances that are known regions of massive star formation. 
Two clumps exhibiting a mass of 3.6 x 10 4 M 0 at a distance 
of 4.5 kpc are in the high mass star forming complex W33 
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Fig. A.2: Relative number distribution of ATLASGAL groups in 
the first quadrant in panel a) and groups in the fourth quadrant 
in panel b) with the maximum velocity difference in a group 
for different input values of the velocity interval to search for 
associated sources in a group, Avj. We use 5 to 15 km s _1 for 
that and a fixed maximum distance between individual sources 
in a group, A d\, of 0.3°. 


(iBieging et al.l 119781) . Our KDA resolution yielding the near 
distance is consistent with trigonometric parall ax m easurements 
resulting in a distance of 2.53 kpc dimmer et al.l 12013 1. The 
complex of sources at 5.5 kpc with high masses between 2 x 10 4 
and 5 x 10 4 M e are located in W51 main. A similar distance of 
5.4 kpc is derived from the trigonometric parallax dSato et al.l 
l2012h (see Appendix IB. 2k A cluster of ATLASGAL clumps at 
a distance around 8 kpc with gas masses higher than 10 4 M e 
are in the W43 complex and W43 south (see Appendix IB. 2k In 
addition, the source with a mass of 5.6 x 10 4 M 0 at 11 kpc is 
G10.47+0.03, an ultracompact HII region dWood & Churchwelj 
1 19891) in the W31 complex (see Appendix lB.2k 


Sources with large distances from the Galactic plane 

To investigate if individual complexes dominate extreme values 
for the distance from the Galactic mid-plane, an illustration of 
its radial distribution is shown for the ATLASGAL sample in the 
first quadrant in the top panel and for sources in the fourth quad¬ 
rant in the lower panel of Fig. Qj] Among the sources, which 
exhibit extremely high distances from the mid-plane indicated 
as red triangles, are the clumps at Ro a i ~ 5.4 kpc and exhibiting 
a height larger than 200 pc below the plane in the upper panel. 
These are located in complex 4 at / «; 6.02° and b ~ -1.23° 
(see Table[3}- Because none of these sources contains an embed¬ 
ded HII region, we cannot analyse the HI absorption and deter¬ 
mine the far kinematic distance from missing HI self-absorption, 
which is more uncertain than using the two methods. The com¬ 
plex is also k nown as NGC 6530, which is a young open clus¬ 
ter, for which iLoktin & Beshenovl (1200 ll) derived a distance be¬ 
tween 560 and 711 pc from Hipparcos trigonometric parallaxes, 
which would locate the region at a smaller height below the 


(a) 



(b) 



(c) 



Fig. A.3: Number distribution of ATLASGAL groups with the 
number of sources per group in a), the velocity dispersion in b), 
and radius in c). 
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Fig. A.4: Correlation plot of the radius and velocity dispersion 
of all groups. The black line indicates the relation by ILarsonl 
(1 198 lh. the green, red, and blue lines show the formula by 
iBallesteros-Paredes et al.l (1201II ) for an FF column density of 
10 21 cm -2 , 10 22 cm' 2 , and 5 x 10 22 cm -2 . The purple line shows 
a power-law fit to the data with an exponent of 0.4 ± 0.02. The 
contour plot is shown in the upper panel with the range of the 
logarithm of the radius and the velocity dispersion divided into 
bins of 0.3. 


disk. We can also use only the HISA method for the clump at 
Real ~ 6.8 kpc with a height of 154 pc below the plane, lo¬ 
cated in complex 262 at l = 49.044° and b ~ -1.078°, for which 
we determine the far kinematic distance. This source is asso¬ 
ciated wit h IRAS 19230+1341 and 6.7-GHz methanol maser 
emission dPestalozzi et al.l[2005l) . Moreover, there is a complex 
of sources at Roai ~ 5 kpc, which are located at a distance of 
158 pc below the mid-plane. One of them at l = 17.608° and 
b = -0.726°, also known as IRAS 18227-1358, is in complex 
75, which consists only of this source. The other two clumps 
at l ~ 18.30° and b * -0.72° are in complex 80, which con¬ 
tains no embedded HII region. We thus locate them at the far 
distance from missing HI self-absorption. Because the displace¬ 
ment from the G alactic p lane of star forming regions is confined 
to ~ 120 pc (iUrquhartet al.l 1201 ll) . the assignement of the far 
distance is unrealistic and the complexes are likely at the near 
distance. In addition, the HISA method is less reliable for the 
far distance than for the near distance. Moreover, the sources 
at a large height above the plane are not associated with 21 cm 


continuum sources either and we also use the HISA method to 
determine the far distance to them. Among those are the clump 
at Real ~ 4.6 kpc with a distance of ~ 185 pc from the disk, 
located in complex 58 at lx 15.029° and b ~ 0.851°, and one at 
Real ~ 5.5 kpc and a height of 178 pc, which is in complex 33 at 
l « 11.92° and b as 0.753°. Approximately 7% of the sources in 
the first and fourth quadrant, which we assign to the far distance, 
have a displacement larger than 120 pc from the plane. We keep 
the far distance for these clumps to use a consistent method to 
resolve the KDA for all ATLASGAL sources. Statistically this 
small fraction does not influence the Galactic distribution. 

Clumps in the fourth quadrant (see lower panel of Fig. IT5l i 
with an extremely high displacement from the plane of 210 
pc and a galactocentric radius of 5.8 kpc are in complex 187 
(7 = 330.044°, b = 1.05 0°). These source s are embedded in the 
IRDC G330.030+1.019 (Peretto & Fulleril200l . which hints at 
the near distance, while we determined the far distance for this 
complex. However, our KDA resolution is uncertain because we 
could use only the HISA method, which is less reliable for the 
far distance solution than for the near distance. Additional com¬ 
plexes, which we locate at the far distance based on this method 
and which are therefore at a large height above the Galactic 
plane are complex 186 (l = 329.718°,/? = 0.804°) with only 
one source at Ro a l ~ 5.8 kpc and a heightof J_60 pc. It is asso¬ 
ciated with the IRDC G329.717+0.785 (IPeretto & Fullerll2009b . 
Other examples are complex 280 (l = 340.088°, b — 0.928°) at 
Real ~ 4.2 kpc and a distance from the plane of 180 p c, em - 
bedded in the IRDC G340.094+0.929 (IPeretto & Fulletillool . 
as well as complex 360 (l = 350.816°, b = 0.514°) atRo a i « 8.4 
kpc and a displacement from the plane of 450 pc, located within 
the IRDC G350.816+0.513 (IPeretto & FulleiOOOl . 


Complexes with extremely large scale heights 


The ATLASGAL sample within a galactocentric radius between 
4.5 and 6.5 kpc shows a flat scale height distribution in Fig. Q/7] 
resulting mainly from a few complexes, which contain several 
sources. Approximately 7% of the sources in the first quadrant 
in the galactocentric radius bins from 4 to 7 kpc are in complex 
261 (/ = 49.2°,/? = -0.7°) (see Appendix IB. 2b with a mean 
Galactic latitude of -0.31°. It is located at a kinematic distance 
of 5.54 kpc and at a mean height of 30.25 pc below the Galactic 
plane. Two additional large complexes in the first quadrant are 
close to the Galactic plane and thus have small heights: Com¬ 
plex 172 (/ = 30.767°,/? = -0.050°) (see Appendix E3 con¬ 
tains ~ 8% of the sources in the first quadrant with a mean lati¬ 
tude of 0.07°, a distance of 8.4 kpc, and a mean height of 10.67 
pc. Approximately 8% of the clumps in the first quadrant are in 
complex 167 (/ = 29.96°,/? = -0.2°) (see Appendix IB. 2 1 ( with a 
mean latitude of -0.04°, a distance of 8 kpc, and a mean height 
of 5.84 pc below the mid-plane. In the fourth quadrant complex 
242 (/ = 336.97,/? = -0.010) consists of ~ 10% of the sources 
within a galactocentric radius between 4 and 7 kpc. It is located 
at a mean latitude of 0.01°, nearby at a kinematic distance of 4.7 
kpc, and thus has a small mean height of 1 pc below the Galactic 
plane. Complex 210 (/ = 333.19°,/? = -0.36°) (see Appendix 
IB .2b has the largest contribution, ~ 12%, to the clumps in the 
fourth quadrant in the bins from 4 to 7 kpc. It exhibits a mean 
latitude of -0.38° and a distance of 3.52 kpc, which results in a 
mean height of 23.37 pc below the mid-plane. 
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Sources with high surface densities 

The effective radius is compared with the gas mass of the AT- 
LASGAL sample in Fig. [23] Some sources are located close to 
the dotted blue line, which indicates a surface density of ~ 1 g 
cm~ 2 . Among those sources with high gas masses concentrated 
to small source radii is G35.2-0.74 with a mass of 423 M e and 
an effective radius of about 0.1 pc, which i s asso ciated with com¬ 
pact HII regions detected at 6 cm bv lDent et alJ ( ll 1984 4. Another 
clump at I = 348.604°, b = -0.912°, which has a mass of 114 
M 0 and a radius of 0.08 pc, also exceeds the surface density of 
1 g creT 2 . It is an UCHIIR , which is also observed in CS (2-1) 
by [Bronfman et al.l (11996h . One source with a higher mass of 
~ 5 x 10 4 M q and a radius of ~ 1.5 pc is located in the W51 
giant molecular cloud complex, associated with W51 IRS1 in 
W51 Main (see Appendix IB.21 ). It emits at 20 pm, 2 cm, and 6 
cm, which shows that it is an exte nded HII region illuminated by 
early OB stars in the centre ( Genzel et all 19821) . Another clump 
with a slightly smaller mass of ~ 3 x 10 4 M Q and effective radius 
of 1.24 pc is in W51 IRS2, which is a compact infrared source 
near W51 IRS1 dGenzel et .al.Hl 9821) . 


Appendix B.2: Identification of complexes with known giant 
molecular clouds 

W43 Main 

Complex 172 has most ATLASGAL sources in the first quad¬ 
rant and consists of 47 clumps . It is located in the well known 
W43 main star fo rming complex (/ = 30.767°, b = -0.050° 
dMotte et al.112003 )'). which contains a giant HII region. The as¬ 
sociation of giant molecular clouds with the W43 stellar cluster 
and the star f ormati on efficiency of the molecular cloud complex 
of ~ 25% dMotte et all2003l) . which might result in a stellar den¬ 
sity of ~ 100 stars pc~ 3 , hint at a sudden ministarburst in W43. 
It can therefore be classified as a Galactic mini-starburst region 
and serve as a template to constrain cloud and gas properties in 
distant starburst galaxies (iMotte et al.ll2003l) . 

W43-south 

We place 30 sources in the first quadrant in comp lex 167 
(/ = 29.96 °,b = -0.2° dNguven Luong et~afl 1 201 in '), which 
is W43-south. It is located at a slightly smaller Galactic lon¬ 
gitude, 29.8° < I < 30.5°, than the W43 main cloud (com- 
plex 172). Investiga t ion o f 13 CO emission toward the two 
(iNguven Luong et al.l 1201 1.1 revealed that they might be con¬ 
nected by low-density gas and that they likely belong to the same 
giant molecular cloud. The distance to the W43 region has al¬ 
ready been derived by several studies: IPratap et al . (1199 m) ob¬ 
tained the near distance to the HII region G29.96-0.02, which 
is located in W43-south, by analysing extinction as well as 
formaldehyde absorption lines. Moreover, Motte et al. (in prep.) 
also derive the near distance to that so urce from large-scale HI 
absorption. lAnderson & BaniaJ (i2009h determine the near dis¬ 
tance to the W43 main cloud (complex 172) from HI self¬ 
absorption, while they locate G29.96-0.02 and another source 
in W43-south at the far distance from their investigation of HI 
self-abs orptio n as well as HI absorption. These methods are also 
used by iRusseil et al.l (1201 ll) . who determine the near distance 
to several sources in the W43 main cloud and note that there is 
a disagreement in distances to G29.96-0.02 given in the liter¬ 
ature. They use the near distance to that source derived from 
the maser parallax method. The W43 complex has also been 


observed by the BeSSeL Survey dReid et al.ll2009a ): From ob¬ 
servations of thre e 12 GHz methanol masers and a 22 GHz 
water maser [Zhang et~afl (i2014l ) derived a parallax distance of 
5.5 ± 0.37 kpc. iKolpak et al. ( 2003b obtained the far distance 
to two sources in W4 3-sou th by analysing HI absorption. Using 
the same method lPandian et al . (200 8) derived the far distance to 
some clumps located in the W43 main cloud. These distance as¬ 
signments agree with our distance resolutions for the W43 main 
cloud (complex 172) and W43-south (complex 167), which we 
place at the far distance. Investigation of HI self-absorption (see 
Sect. 15.2b toward the two complexes hints at the near distance, 
while HI absorption (see Sect. 15.4b toward four sources in the 
W43 main cloud and two clumps in W43-south reveals the far 
distance. Combining the two methods as described in Sect. 15.51 
gives the far distance of 8.4 + 0.82 kpc to complex 172 and 
7.9 + 0.64 kpc to complex 167. 


M17 


We locate 22 sources in complex 50, which is M17 (the Ome ga 
Nebu la, the Swan Nebula) (I = 15.10°, b = -0.75° ( Wu et al.l 
(2009)). It is one of the giant molecular_clouds_with the largest 
sizes and masses in the Milky Way ( Ladai 11976h . A study of 
dust em issio n _at near-infrared, far-infrared, and submm wave¬ 
lengths dChini et al.lll987f) showed that M17 is associated with 
one of the most lumino us compact HII regions in the Galaxy. 
lLada & Chaissonl (i 1 9751 ) observed 6-cm H 2 CO aborption lines 
toward Ml7, which revealed a lane of dust obscuration west of 
the HII region. This part of the molecular cloud contains a small 
region of strong CO emission located south-west of the contin¬ 
uum peak at 6 cm (Ml7 SW) and is also detected in other molec¬ 
ular lines such as HCN and H 2 CO at 140 GHz (lLada et al.ll 1974f) 
probi ng high cloud densities of > 10 5 cm -3 ( Thaddeus et aD 


11971 ). Varying H 2 CO linewidths ( Lada & Chaissonlll975l) and 
CO measurementsjn M17 SW, which give a lower mass limit 
of 6000 M q ( Lada et al . 1974 :). hint at a molecular cloud col¬ 
lapse. The 12 G Hz CH;OH ma ser line emitted in M17 SW 
was measured by IXu et all (1201 ll) to determine the trigonomet¬ 
ric parallax. They obtain a distance of 2 + 0.13 kpc and locate 
G15.03-0.68 in the Sagittarius spiral arm. We also derive a near 
distance of 2.4 + 0.61 kpc from HI self-absorption and HI ab¬ 
sorption toward two sources in the complex. 


W51 Main 

Complex 261 with 22 ATLASGAL cl umps is part of the W51 
molecular cloud (/ = 49.2°, b — -0.7° (lGreenil2009l) ). known as 
W51 Main. The whole W51 complex is one of the most mas¬ 
sive (1.2 x 10 6 M g ) a nd largest (Al x A b = 83 X 114 pc) gi¬ 
ant molecular clouds ( {Carpenter & Sanderslll998l) . In addition, 
W51 is one of the most luminous (10 7 L 0 ) s tar forming com¬ 
plexes and harbours many O stars (iBiegine 1975 ). which in¬ 
dicate that the region is currently forming an OB association. 
Moreover, iKumar et al.l 0 2004 ) discovered embedded star clus¬ 
ters from their near-infrared observations, which are associated 
with UCHIIRs. They obtain a high star formation efficiency of 
~ 10% f rom an estimate of the W51 molecular cloud mass. 
i Sato et al.l (120121) observed the 22 GHz H 2 O maser line in W51 
Main to measure the trigonometric parallax, which results in a 
distance of ~ 5.41 kpc. The 21 cm continuum data toward four 
ATLASGAL sources show absorption lines at velocities between 
65 and 73 km s” 1 , which are consistent with the higher velocity 
range of some of the absorption spectra in the W51 molecular 
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cloud by iKool dl997l) . Previous studies placed the complex in 
the S agittarius arm at the tangent point dKooll 1 9971 : lEtoka et al. 
1201 2l) . which agrees with our distance assignment at 5.5 ± 1.67 
kpc. 


W31 


We place 21 sources in comp le x 22, which is W31 (/ = 
10.3262°, b = -0.1432° (iGoss & Shaved fl970h ). It is a mas¬ 
sive star forming region, which was resolved into two com¬ 
pon ents (G10.2-0.3 and G 10.3-0.1) by IGoss & Shaved (fl970h 
and lShaver & Gossl (Il970t> using observations at 5 GHz and 408 
MHz. Complex 22 also includes two UCHIIRs (G10.15-0.34 
and G10.30-0.15. IWood & Churchwelli IT989I) located in the 
larger scale surrounding of the two subregions. iDownes et al] 
( 1980 ) compared velocities of the H 11 0a hydrogen recombina¬ 
tion lines of clumps in W31 with HiCO absorption line veloci¬ 
ties to resolve the KDA, which results in the near kinematic dis¬ 
tance of 6 kpc. lBlum et al. d200ll) measu re a spectrophotometric 
distance of 3.4 kpc, while Wilsonl d 19721) obtains H 2 CO absorp¬ 
tion lines at velocities, which are higher than the hydrogen re¬ 
combination line velocity, and therefore assig ns the complex to 
the far distance. iGeorgelin & Georgelinld 19761) also measure ab¬ 
sorption at velocities exceeding the velocity of the HII regions, 
which hints at the far distance of 18.7 kpc. In addition, analysis 
of H 2 CO, OH, and HI absorption lines by ICorbel et al.l d 19971) 
indicates the far kinematic distance of 14.5 kpc. Examination of 
HI self-absorption toward the complex yields the near distance. 
In addition, HI absorption toward two clumps (G10.15-0.34 and 
G10.32-0.16) is at velocities, which are larger than the source 
velocity, but much smaller than the tangent point velocity and 
therefore hints at the near distance of 2 ± 0.93 kpc. However, 
this analysis cannot reject the far distance because absorption 
at velocities larger than ~ 50 km s“ ' is unlikely given the 
lack of HI and CO ga s at Galactic longitudes between 5° and 


25° ( s ee e.g. Fig. 3 in Dame et aD 1 19871 ICorbel & Eikenberrvl 


l2004 . ISanna et al.l (hOllT observed 22.2 GHz H 2 0 maser emis¬ 


sion toward G10.62-00.38 t o measu re the trigonometric parallax 
within the BeSSeL Survey dReid et al.l2009al) which resulted in¬ 
deed in a distance of 5 kpc and located W31 on the near edge of 
the hole in the gas distribution within a galactocentric radius of 
~ 3 kpc. 


N49 

Complex 156 consisting of 17 ATLASGAL clumps is asso¬ 
ciated with the HII region N49 (l = 28.827°, b = -0.229° 
dChurchwell et alJl2006l) ). It is one of the bubbles, which are de- 
tected by the investigation of GLIMPS E dBeniamin et al Jl2()03l) 
images from IChurchwell et al.l d2006l) . This region consists of 
a layer of dense neutral material, which accumulates as a shell 
around a radio HII region with an approximately spherical shape. 
20 cm free-free emission probes the ionized gas of N49, the 
thermal continuum at 24 fim traces the hot dust within the HII 
region, which is_ surrounded by an outer layer of 8 fim emis¬ 
sion d Watson et al .1120081) . A hole is revealed in the centre of 
the 24 /urn and 20 cm emission, where the exciting star is lo¬ 
cated. This mig ht be an 05 V or 08111 star, although the age es¬ 
timate by 1 Watson et al.l d2008i) favours the 05 V star with > 10 5 
yr. Analysis of the 870 /urn emission yields four condensations 
within the shell around the central HII region and a mass of 
4200 M 0 for the whole shell dDeharveng et aljf2010t) . NH 3 ve¬ 
locities dWienen et al.ll 1201 2 l) are observed toward three con¬ 


densations and the 13 CO velocity dJackson et alJl2006b is mea¬ 
sured toward the fourth condensation indicating their assoca- 
tion with the ionized gas. Two massive young stellar objects 
can be revealed at 24 pm , which are in an early evolutionary 
phase (IWatson et al.il2008h . as well as an UCHIIR l ocat ed in the 
surro u nding of the most massive condensation dDeharveng et al.l 
l2010h . lAnderson & Bani a (2009) derived the near distance to the 
region from HI self-absorption and HI absorption. Because none 
of the observed ATLASGAL sources in N49 satisfied our criteria 
for radio continuum absorption, we could only look for HI self¬ 
absorption and we also determined the near distance of 5.7+0.41 
kpc. 


M16 

We locate ten sources in complex 70, which is M 16 (the Ea¬ 
gle Nebula) (/ = 16.94°, b = 0.76° dZeilik & LadallT97l ). The 
HII region M16 belongs to a larger star forming complex, which 
is mo stly still embedded within the molecular cloud dGoudisI 
119761) . Images taken by the Hubble Space Telescope reveal struc¬ 
tures in the molecular cloud containing dense globules, which 
are photoevaporated more slowly than their lower density envi¬ 
ronment, forming ’’molecular pillars”. The heating and ioniza¬ 
tion of M16 results from the young open star cluster NGC6611. 
The Eagle Nebula is associated with recent and ongoing high 
mass star formation. Its stellar population has been analysed by 
iHillenbrand et alJ d!993l) . They estimated an age of (2 ± 1) x 10 6 
yr and locate the region at a distance of 2 kpc from spectro¬ 
scopic parallaxes. To resolve the KDA we could only analyse 
HIS A because complex 70 contains no 21 cm continuum source. 
We obtain a near dista nce of 2.2 ± 0,56 kpc, wh ich agrees with 
the distance derived by IHillenbrand et akl (119931) . 


G333 

Complex 210 consists of most ATLASGAL sources in the fourth 
quadrant containin g 58 clumps (l = 333.19°, b = -0.36° 
dBains et al.l 120061) '). It is known as the G333 giant molecular 
cloud complex or RCW106 region with several molecular clouds 
forming high mass stars as well as HII regions. We determine the 
near kinematic distance to complex 210 of 3.5+0.35 kpc from HI 
self-absorption and HI absorption toward several sources. Our 
estimate is con siste nt with the kinematic distance of 3.6 kpc cal¬ 
culated by iLockm anl (i 19791) as well as the range from 3.3 to 3.8 
kpc given bv lUrauhart et al.l ( 1201 2 l) an d with the spectroph oto- 
metric distance of 3.96 kpc obtained bv lMoises et al. (1201 li) . 

G305 

From our ATLASGAL sample in the fourth quadrant 41 sources 
are located in comp lex 29 (l = 305.361°, b = 0.056° 
dClark & Ported 1 2004 )). which is the rich high mass star form¬ 
ing complex G305. Some luminous HII regions are embed¬ 
ded within the _molecular cloud complex, the analysis of MSX 
data by IClark & Ported (120041) yields strong ra dio sourc e s with 
up t o 31 07V s ta rs ion izing an HII region. iRussei 1 d2003t) 
and IMoises et al.l (1201 ll) divided the region into two parts, 
G305.2+0.0 and G305.2+0.2, and derived distances to them. 
While lRusseill (120031) obtained a distance o f 3.5 kpc for b oth, the 
spectrophotometric distance analysed by IMoises et alJ (1201 ll) 
is inconclusive for G30 5.2+0 .0, but yields the same result 
as assigned by iRusseill (120031) for G305.2+0.2. Our HI self¬ 
absorption method also indicates the near distance to complex 
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G305. However, the analysis of HI absorption toward one source 
of the region, G305.19+0.03 (see Fig. IB. II ). results in the far dis¬ 
tance. The composite of the two methods reveals the far distance 
of 6.6 ± 1.24 kpc to the complex. 

NGC 6334 

Complex 356 with 28 ATLASGAL clumps is known as the 
NGC 6334 molecular c loud complex (l = 351.317°, b = 0.661° 
dRodriguez et al Jfl 982h ) with embedded HII re gions identified 
at 6 cm as compact to extended radio sources. iNeckell (Il978h 
derived the distance to NGC 6334 from optical observations of 
young stars within this complex of 1.74+0.31 kpc and place 
it in the Sagittarius arm. They also give a kinematic distance 
of 0.7+1.9 kpc, which corresponds to the interarm region. The 
two values are within the uncertainty of the near kinematic dis¬ 
tance of 0.8 ± 1.03 kpc, which we obtained. The large distance 
uncertainty originates from the unreliability of kinematic dis¬ 
tances near a longitude of 0°. Our result also agrees with the 
spectr ophotometric distance of 1.82 kpc given bv lMoises et al.l 
(1201 ll) and with the distance of 1.35 k pc derived r ecentl y from 
trigonometric parallax measurement bv lReid et al.l (2014 !). 


G305.19 + 0.03 



Fig. B.l: HI spectrum toward G305.19+0.03, an HII region in 
the high mass star forming complex G305 extracted from the 
SGPS. For explanation of the different lines see Fig. [3] The HI 
absorption method yields the far distance. 

Appendix C: Comparison of kinematic and 
trigonometric parallax distances 


G338.4+0.12 


Complex 255 contains 24 ATLASGAL sources, which con¬ 
nects t hree small complex e s. The complex is a giant HII region, 
where ICaswell & Havnesl ( 19871) measured h y droge n recombi¬ 
nation lines near 5 GHz. ICaswell & Havnesl (Il987l) derived a 
vast range of recombination line velocities over 120 km s _1 in¬ 
dicating several features, which lie on the same line of sight to 
the observer at different distances. We obtain a narrower velocity 
range from -30 to -42 km s 1 of most sources and a few clumps 
exhibit - 46 to -56Jem s 1 in our smaller complex at a longitude 
of 338°. iGeorgelin et al . ( i 19961) also determined velocities and 
distances to two radio sources in the complex from observations 
of Ha emission. They give kinematic distances of 2.2 kpc for 
338.398+0.164 and 3.2 kpc for 338.450+0.061, which lie within 
the errors of our distance estimates. However, they should as¬ 
sign the far distance^ taking the OH absorption line into account. 
While lRusseill ( 2003 ) also places the complex at the far distance, 
our KDA resolution reveals the near distance of 3.06 kpc. The 
HI intensity map shows absorption, which hints at the near kine¬ 
matic distance. In addition, the HI absorption method results in 
the near kinematic distance to one HII region and the far dis¬ 
tance to four 21 cm continuum sources. The combination of the 
two techniques gives the near kinematic distance of 3.1 ± 0.41 
kpc to the complex. 


Gum 50 

We locate 23 sources in the fourth quadrant in complex 163, 
which is known as Gum 50 (l = 328.573°, b = -0.531° 
dBeichman et al JlT 988l) ). Observations of hydrogen recombina- 
tion lines near 5 GHz by ICaswell & Havnesl (1 19871) reveal HII 
regions located in the complex such as IRAS 15539-5353 and 
G328.57-0.53. Moreover, a few sources in this high mass star 
forming region, e.g. IRAS 15525-5407, IR AS 15539-5353 
and G 328.34-0.53, are associated with IRDCs dPeretto & Fullerl 
120091) . which indicates the near distance. This is in agreement 
with our HI self-absorption and HI absorption analysis resulting 
in the near kinematic distance of 3 ± 0.39 kpc. 


Recently, more than a hundred distances to high mass star form¬ 
ing regions have been determined via measurements of trigono¬ 
metric parallaxes of class II CH 3 OH and H 2 O masers with Very 
Long Baseline Interferometry (VLBI). The data were taken with 
the VLBI Exploration of Radio Astronomy (VERA) arra\0 and 
with the NR AO Very Long Baseline Array in the cours e of 
Bar and Spiral Structure Legacy (BeSSeL) surve\0 dReid et al.l 
120141) . Parallaxes provide gold standard distance measurements 
since they result from direct measurements and basic trigonom¬ 
etry and do not require assumptions such as a rotation curve. 

Of the 103 sources for which lReid et al. ( 2014 ) present par¬ 
allax distances, 48 are covered by ATLASGAL. For total of 27 
of these we have determined kinematic distances. Our LSR ve¬ 
locities, determined from NH 3 spectra agree wit h ma ser veloc¬ 
ities within the uncertainties assigned bv iReid et al.l ( 2014t) . in 
15 cases within 2 km s _1 and in all cases within 6 km s -1 . In 
Fig- EH we compare the parallax distances with our kinematic 
distances. 

For a third of the 27 sources the kinematic and the parallax 
distance agree within the errors and all but two within a factor of 
two. For the latter two, G12.88+0.48 and G16.58-0.05, choos¬ 
ing the far kinematic distances, 12.9 and 11.7 kpc, respectively, 
obviously was a wrong decision, which for the former is rein¬ 
forced by its high latitude. Both are found in complexes that do 
not contain radio continuum sources and thus have no HI absorp¬ 
tion. They were placed at the far distance because they also did 
not exhibit HI self-absorption. 

Part of the scatter in Fig. IC.ll mav be caused by peculiar ve¬ 
locities, which of course are not accounted for in the rotation 
curve model used to determine the kinematic distances. A dra¬ 
matic example for this is the case of W30H for which the paral¬ 
lax distanc e is al most a factor of two smaller than the kinematic 
distance (IXu et al.1120061) . 

Some of the sources with measured parallax distance and re¬ 
solved kinematic distance ambiguity are described in more de¬ 
tail. 


2 http://veraserver.mtk.nao.ac.jp 

3 http://bessel.vlbi-astrometry.org 
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Table B.l: Associations of largest ATLASGAL complexes with known giant molecular clouds (see Sect. 17. II) . 


Sample 

Complex number 

Number of clumps 

Giant molecular cloud 

Distance C0mp iex 

(kpc) 

log 10 (Sum mass clumps /M 0 ) 

in first quadrant 

172 

47 

W43 Main 

8.4 

5.8 


167 

30 

W43-south 

7.92 

5.54 


50 

22 

M17 

2.36 

5.14 


261 

22 

W51 Main 

5.54 

5.68 


22 

21 

W31 

1.97 

4.67 


156 

17 

N49 

5.71 

4.84 


70 

10 

M16 

2.22 

4.32 

in fourth quadrant 

210 

58 

G333 

3.52 

5.49 


29 

41 

G305 

6.58 

5.52 


356 

28 

NGC6334 

0.83 

4.09 


255 

24 

G338.4+0.12 

3.06 

4.75 


163 

23 

Gum 50 

2.96 

4.61 



0 


5 


10 


Trigonometric Parallax Distance (kpc) 


Fig. C.l: Kinematic distances determined in our study vs. 
trigonometric parallax distances from the BeSSeL survey and 
VERA results. The dotted fiducial line marks equality. 


G12.89+0.49 


HIS A method is less reliable for the far distance than for the near 
distance. 


G12.91-0.26 and G12.68-0.18 

These two sources are located in the W33 molecular cloud 
comp lex. It is a strong radio continuum source detected by 
IWesterhoutt (i 19581) and consists of various evolutionary phases 
of high mass star formation from quiescent IRDCs to HII re¬ 
gions. Radio recombination line observations revealed two ve¬ 
locity components in different parts of W33, one at ~ 36 km s 1 
and one at ~ 58 km s~ ! dBieging et al.l!978l) . We also obtain sev¬ 
eral velocities from NH 3 observations in this region and divide 
it into different complexes. We locate G12.91-0.26 in complex 
38 with a mean radial velocity of 43 km s and a velocity dis¬ 
persion of 6.81 km s -1 and G12.68-0.18 in complex 41 with a 
mean velocity of 55.56 km s~' and a velocity dispersion of 0.57 
km s -1 . llmmer et al.l (1201 3l) used a water maser exhibiting ve¬ 
locities between 36 to 38 km s _1 to derive the parallax of their 
observation at l = 12.909°, b = -0.261°, which corresponds to 
a distance of 2.53 kpc. This is consistent with our investigation 
of HI self-absorption and HI absorption towards two sources in 
comp lex 38, which reveal the near distance of 4.2 ± 0.39 kpc. 
llmmer et al.l (1201 3l) also fit the position of a water maser to ob¬ 
tain the parallax of G12.68-0.18, which results in a distance of 
2.4 kpc. This is in agreement with our KDA resolution of com¬ 
plex 41, which gives the near kinematic distance of 4.9 + 0.3 kpc 
using only the HISA method. 


G12.89+0.49 is in comp lex 44. This source is a high mass 
protostellar object dBeuther et al. 2002 ). which is associated 
with IRAS 18089-1732 (ISridharan et aD 120021). It is well 
known to harbour s trong HoO (lAng lada et alT 19961). C H;OH. 
(ICaswell et a1 .11 1 9951) and OH mas ers rtC ohen et all 19911) . Previ¬ 
ous studies ( Sridharan et al JI2002L iBeuther et al. 20021) give the 
near and far kinematic distance, but did not solve the KDA. Us¬ 
ing 12.2 GHz methanol masers IXu et al.l (1201 ll) derived a dis¬ 
tance of 2.34 kpc from the trigonometric parallax, llmmer et al.1 
( 2013 i) also measured a parallax distance from 22.2 GHz wa- 
te r maser obse rvations, which is in agreement with the result of 
IXu et al.l! ( 201 ll) . and locate the high mass star forming region 
in the Scutum spiral arm. Because there is no 21 cm continuum 
source in complex 44, we could only analyse HISA, which re¬ 
veals the far kinematic distance of 13 + 0.46 kpc. However, the 


G23.01 -0.41 

G23.01-0.41 is located in complex 114. This high mass star 
forming region contains a 12 GHz methanol maser, which co¬ 
incides with a H 2 CO maser and NH 3 (3,3) peak emission. A 
fit to the positions of the two masers gives one parallax, which 
corresponds to a distance of 4.59 kpc dBrunthaler et al.l[2009h . 
Different distances to th i s sour c e are given by previous stud¬ 
ies so far : ICodella et al.l (1 19971) . iForster & CaswelJ (1 19991) and 
ICaswell & Havnesl (1 19831) used the far kinematic distance, al¬ 
though the dis_tance ambiguity was not resolved. In contrast, 
lHariu et all (119981) locate the region at the near kinematic dis¬ 
tance. This is consistent with our KDA resolution resulting in 
4.9 ± 0.29 kpc from HI self-absorption. 
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G23.46-0.18 

G23.46-0.18 belongs to complex 117. This is a star forming 
region, where 6.7 GHz methanol maser emission observed by 
I Walsh et~aP(ll998l) hints at massive young stellar objects. Previ¬ 
ous studies revealed different velocity c omp onents of the molec¬ 
ular cloud harbouring G23.44-0.18: lOhishi et al.1 d2012h ob¬ 
tained three component s at 101, 1 03, and 104 km s 1 from their 
H 13 CO + spectra, while iLockman ( 198 9) derived the radio re¬ 
combination line velocity of 103 km s for an HII region, which 
is likely associated with the molecular cloud. iRathborne et al.l 
(l2009h give two 13 CO velocities at 101 km s 1 and ~ 104 km s 1 , 
for which iRoman-Duval et alJ (120091) determined kinematic dis¬ 
tances of 6.43 kpc to the component at 101 km s -1 and of 6.65 
kpc to the 104 km s -1 com ponent from HIS A and HI absorp- 
tion. iBrunthaler et~ahl d2009h fitted the positions of four 12 GHz 
methanol masers to measure the parallax of their observation 
at 1 = 23.4398°, = -0.1822°. This yields a distance of 5.88 

kpc and places the region near the end of the galactic bar in 
the Norma arm. We measured a velocity of 98.4 km s 1 toward 
G23.46-0.18, which is clo se to the vel ocity of the methanol 
maser at 97.6 km s 1 ( iBrunthaler et al . 200 9:). We derive the 
near kinematic distance of 5.9 ± 0.3 kpc using the HIS A method, 
which is in agreement with the distance from the trigonometric 
parallax and also similar to the kinematic distance of 6.43 kpc of 
IRoman-Duval et aTl d2009l) . 

G23.96-0.il 

G23.96-0.ll is in complex 121, which is associated 
with G23.657-0.127 observed with the BeSSeL survey. 
G23.6 57-0.127 is a massive protostar, where IBartkiewicz et al.l 
(120051) detected 6.7 GHz methanol maser emission between 
77 and ~ 88 km s _I . Analysis of an EVN image showed 
that the met hanol masers are arranged in an approximately 
circular ring (IBartkiewicz et al][2005l) . Bartkiewicz et al.l (120081) 
observed 12.2 GHz methanol masers with the VLBA towards 
G23.657-0.127, which have all 6.7 GHz counterparts and are 
also distributed spherically symmetric. Fitting 19 positions of 
12.2 GHz methanol masers, which have approximately the same 
velocities as the 6.7 GHz methanol masers, reveals_a distance 
of 3.1 9 kpc from the trigonometric parallax (IBartkiewicz et alJ 
l2008l l. We obtain a similar velocity range of complex 121 with 
a mean velocity of 78 km s 1 and a velocity dispersion of 3.7 
km s _1 . The HIS A method and HI absorption towards one 
source in complex 121 results in the near kinematic distance 
of 4.9 ± 0 .31 kpc , which is similar to the distance derived by 
IBartkiewicz et al.l(l2008l) . 

G35.19-0.74 

G35.19-0.74 is located in complex 214 and associated with 
G35.2-0.7, which is observed by the BeSSe L sur vey. The com¬ 
plex is a massive star forming region, where iGibb et aTi (120031) 
observed CO, SiO, and radio emission indicating several out¬ 
flows from G35.2-0.7. In addition, VLA observations at 3.5 and 
6 cm showed fragmentation of the radio continuum into mul- 
tiple YSOs embe d ded within den s e clum ps along the outflow 
(iGibb et alJ 12003 ). I Solomon et al.l dl987l) derived the near dis¬ 
tance of 3.3 kpc because the far distance would locate the source 
at a larger height than 150 pc from the plane. The positions of 
two 12 GHz m asers w ith radial velocities of 27.9 and 27.5 km s 1 
were fitted by IZhang et~akl d2009l) to obtain the trigonometric 
parallax, which leads to a distance of 2.19 kpc. The mean ve¬ 


locity of complex 214 is 34.1 km s -1 with a velocity dispersion 
of 2.02 km s _1 , which agrees with the velocities of the masers. 
Because complex 214 contains no 21 cm continuum source, we 
use only the HISA method, which reveals the near kinematic 
distance of 2.3 ± 0.43 kpc. Our KDA resolution is therefore con- 
sistent with the di stances given by ISolomon et al.l d!987l) and 
IZhang et~akl d2009l) . 
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